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Abstract 


This  report  presents  parts  of  the  scientific  investiga¬ 
tions  carried  out  under-  a  task  program  related  to  the  deter¬ 
mination  and  assessment  of  surface  dynamic  soil  properties 
by  laboratory  instrumentation.  It  first  reviews  the  histor¬ 
ical  development  of  the  hardness-penetrometer  concept  and 
presents  a  pilot  model  for  a  new  aerial  penetrometer.  Basic 
theoretical  considerations  are  then  advanced  regarding  the 
stress-strain  relationships  in  axial  symmetry  and  the  dynamic 
loads  associated  with  a  drop-impact  penetrometer.  The  exist¬ 
ing  theory  of  elasticity  solutions  are  then  reviewed  and  the 
solutions  for  a  sphere  impacting  on  an  ideally  plastic  and 
on  an  elasto-plastic  semi-infinite  solid  are  presented.  The 
Hertz  solution  of  an  elastic  sphere  impacting  on  an  elastic 
half-space  is  modified  by  using  the  elastic-viscoelastic 
analogy  to  find  a  solution  for  a  sphere  dropping  onto  a  four 
parameter  viscoelastic  half-space. 

The  problem  of  energy  transfer  during  impact  by  vibra¬ 
tions  is  discussed  shortly.  The  data  for  fifteen  test  series 
carried  out  with  the  pilot  model  penetrometer  on  natural  soils 
and  on  some  artificial  materials  are  presented  and  evaluated. 
It  is  concluded  that  the  principle  of  the  drop  impact- 
penetrometer  is  applicable  in  experiments  for  assessing  the 
mechanical  properties  of  a  natural  soil  surface.  The  proposed 


modification  of  a  prototype  penetrometer  should  also  increase 
its  capability  and  reliability  as  a  remote  sensor  for  opera¬ 
tional  applications.  The  instrument,  as  envisaged,  will  not 
only  permit  determination  of  mechanical  surface  properties 
on  terrestrial  soils  but  should  also  be  extremely  useful  for 
the  exploration,  by  remote  techniques,  of  lunar  and  other 


extraterrestrial  surfaces. 
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I •  Introduction  and  Definition  of  Problem 

The  determination  of  soil  properties  is  one  of  the  most 
important  tasks  of  modern  soil  mechanics.  For  its  successful 
solution  one  needs,  first  of  all,  a  clear  and  rigorous  defini¬ 
tion  and  a  thorough  understanding  of  these  properties.  Next, 
it  will  be  necessary  to  develop  suitable  methods  and  instru¬ 
ments  for  determining  these  properties  and,  finally,  one  will 
have  to  evaluate,  interpret  and  apply  the  data  that  are 
obtained . 

Historical  progress  in  this  field  followed  a  more  or  less 
logical  development.  It  began  by  measuring  properties  of  dis¬ 
turbed  soils  in  the  laboratory  (e.g.  Atterberg  consistency 
limits,  friction  angle  in  the  shear  box).  But  one  soon  recog¬ 
nized  that  certain  soil  properties  were  strongly  influenced 
and  adumbrated  by  the  disturbance  that  is  unavoidable  in  the 
normal  soil  sampling  process.  For  this  reason  the  measurement 
of  properties  in  the  laboratory  on  "undisturbed  samples"  found 
wider  and  wider  acceptance  (e.g.  the  consolidation  test  and 
the  triaxial  test).  It  is  obvious  that  the  term  "undisturbed" 
must  be  understood  relatively  and  merely  designates  generally 
accepted  standard  methods  by  which  soils  are  being  sampled  to 
achieve  reasonably  undisturbed  samples  without  undue  expense. 
The  actual  method  used  for  sampling  and  its  applicability  for 
the  soil  being  sampled  becomes  very  important. 


Even  these  methods,  however,  do  not  yet  constitute  an 
ideal  solution  because  first,  a  truly  undisturbed  sample  can¬ 
not  be  obtained  -  even  with  the  greatest  care,  and  secondly, 
the  transport  of  the  sample  to  the  laboratory,  the  assembly 
of  the  sample  in  the  testing  apparatus,  the  test  itself  as 
well  as  its  evaluation  always  require  a  large  amount  of  time 
and  labor.  Moreover,  according  to  Burmister  (1957)  the  soil 
properties  depend  also  on  the  soil  environment,  particularly 
the  system  of  stresses  that  exists  "in  situ"  in  the  soil  and 
this  author  is  in  complete  agreement  with  this  hypothesis. 

For  this  reason  it  would  be  advantageous  to  develop 
methods  that  permit  the  measurement  of  soil  properties  "in 
situ,"  i.e.,  in  place.  The  advantages  of  such  methods  could 
be  listed  as  follows: 

1)  There  is  no  need  for  sampling  and,  therefore,  any  dis¬ 
turbance  during  the  sampling  process,  the  transport 

of  the  sample  and  its  assembly  in  the  apparatus  is 
eliminated . 

2)  The  time  lag  between  sampling  and  evaluation  of  the 
test  results  is  eliminated. 

3)  The  measured  properties  are  available  immediately  on 
location. 

For  these  reasons  developments  of  new  measuring  techniques 
have  been  directed  increasingly  towards  instruments  and  methods 
that  measure  soil  properties  "in  situ."  Among  these  instru¬ 
ments  and  methods  are: 


a)  The  standard  penetration  test 

b)  The  shear  vane 

c)  The  Menard-pressuremeter 

d)  The  static  penetrometer  (e.g.  C3R) 

e)  The  drop  penetrometer  (e.g.  aerial  penetrometer) 

The  methods  a),  b)  and  c)  are  mainly  concerned  to  meas¬ 
ure  the  shearing  strength  of  a  soil  in  a  bore  hole  whereas 
methods  d)  and  e)  have  the  purpose  to  determine  the  properties 
of  the  soil  surface.  The  mechanical  properties  of  the  soil 
surface  are  of  importance  in  the  design  of  highway  or  airport 
pavements,  or  when  it  is  necessary  to  judge  the  trafficabil- 
ity  of  a  surface  by  a  tractor,  by  construction  equipment  or 

by  a  military  vehicle. 

The  latter  problem  is  a  tactical-logistical  one  of  great 
military  significance  which  normally  is  made  more  difficult 
still  because  often  the  soil  surface  to  be  judged  is  not 
accessible.  From  a  purely  scientific  point  of  view  this  prob¬ 
lem  is  of  great  interest  also.  For  example,  the  design  of 
the  lunar  exploratory  module  (LEM)  being  developed  at  the 
moment  by  the  Grumman  Aircraft  Company  requires  a  reliable 
indication  of  the  surface  properties  of  the  moon. 

In  the  fall  of  1962  the  author  initiated  a  research  pro¬ 
ject  sponsored  by  the  Terrestrial  Sciences  Laboratory,  Air 
Force  Cambridge  Research  Laboratories  (OAR),  Hansccm  Field, 
Bedford,  Massachusetts.  A  part  of  this  project  was  to 
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investigate  the  principles  which  would  permit  the  measurement 
of  mechanical  properties  of  a  soil  surface  by  remote  methods, 
i.e.  without  requiring  the  physical  presence  of  an  operator 
at  the  surface.  The  aim  was  to  determine  the  properties  of  a 
natural  soil  surface  that  might  control  modern  aerospace  oper¬ 
ations  on  natural  terrain  such  as  landing  and  take-off  of  air¬ 
craft  and  recovery  of  space  vehicles. 

The  problem  then  was  to  develop  a  method  for  measuring 
those  properties  of  a  soil  surface  which  govern  the  opera¬ 
tional  possibilities  of  natural  terrain  sites  for  modern  mili¬ 
tary  aircraft  and  other  aerospace  vehicles.  The  problem  was 
not  .entirely  new. 
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During  1950-1953  an  aerial  penetrometer  was  developed  at 
New  York  University  (Fig.  1;  see  also  Molineaux,  1955*  and 
Anonymous,  1953)  that  could  be  dropped  from  an  airplane  and, 
upon  impact,  would  release  signals  in  the  form  of  flares, 
colored  vanes  or  flags  that  would  indicate  the  severity  of 
the  impact.  The  signals  could  be  observed  from  the  plane  and 
would  indicate  the  hardness  of  the  impacted  surface. 

In  its  practical  application,  this  aerial  penetrometer 
had  several  flaws.  One  of  the  principal  difficulties  was  that 
it  required  an  impact  with  zero  horizontal  velocity  component 
together  with  a  perfect  vertical  alignment  of  the  penetrometer 
axis  . 

In  the  normal  operational  case  this  was  to  be  achieved 
by  flying  the  airplane  at  a  constant  speed  of  180  knots  and 
ejecting  the  penetrometer  rearwards  and  slightly  downwards 
with  a  speed  such  that  the  net  difference  between  the  two  hor¬ 
izontal  velocity  components  would  be  zero.  Because  of  wind 
and  small  variations  in  the  forward  speed  of  an  airplane  it 
is  most  difficult  to  achieve  exactly  this  condition  of  =  0, 
Even  more  important,  if  the  longitudinal  axis  of  this  pene¬ 
trometer  is  just  slightly  off  the  vertical  at  the  instant  of 
impact,  part  of  its  energy  will  be  converted  into  a  moment 
which  tends  to  throw  the  penetrometer  onto  the  surface  broad¬ 
side  and  thus  a  good  part  of  the  impact  resistance  will  not 
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be  transmitted  through  the  point  of  the  penetrometer.  It  was 
found  rather  difficult  to  stabilize  the  penetrometer  aerody- 
namically  in  such  a  fashion  that  it  would  always  hit  verti¬ 
cally.  Furthermore,  this  instrument,  as  developed,  cannot 
make  an  absolute  measurement.  Rather,  before  being  dropped  an 
estimate  must  be  made  with  respect  to  the  expected  impact 
resistance  and  the  observed  signal  then  merely  indicates 
whether  the  actual  resistance  encountered  was  larger  or 
smaller  than  the  one  previously  estimated. 

Principle  of  Penetrometer  Methods 

The  development  of  the  aerial  penetrometer  just  discussed 
was  not  at  all  without  precedent  but  followed  in  principle  the 
"cone  penetrometer"  developed  by  the  Waterways  Experiment 
Station  of  the  Corps  of  Engineers,  U.S.  Army.,  in  Vicksburg, 
Mississippi  (Tig.  2)(WES  1957)  which  is  used  to  determine  the 
trafficability  of  a  natural  soil  surface. 

This  instrument  is  a  rod  with  a  standard  conical  tip  at 
the  lower  end.  This  tip  is  pushed  vertically  into  the  soil 
surface  by  a  force  on  the  handle  at  the  upper  end  and  this 
force  is  measured  bj'  a  dial  through  a  proving  ring. 

The  instrument  and  its  operation  is  described  in  a  Water¬ 
ways  Experiment  Station  manual  and  is  discussed  by  Hechtl 
(1964)  in  the  context  of  trafficability  evaluation.  Although 
the  Vicksburg  cone  penetrometer  is  an  original  instrumentation 
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development,  it  is  based  on  the  same  principle  as  other  pene¬ 
trometers.  In  the  final  analysis  all  penetrometer  methods 
have  their  origin  in  the  principles  of  hardness  measurement 
of  materials  which  can,  according  to  Williams  (1942),  be 
traced  back  to  an  article  by  Reaumur  in  1722.  The  more  modern 
investigations,  however,  most  probably  were  initiated  as  a 
consequence  of  the  solution  of  The  contact  stresses  of  two 
elastic  spheres  by  Hertz. 

The  idea  to  use  the  results  of  this  solution  for  the 
determination  of  material  properties  was  obvious  and  led  to 
the  various  methods  of  measuring  the  surface  hardness  of  a 
material. 

Hardness  Measurements 

In  measuring  the  well-known  Brinell  hardness,  e.g.,  a 
steel  sphere  is  pressed  by  a  force  P  into  the  surface  to  be 
examined  and  the  diameter  of  the  resulting  permanent  deforma¬ 
tion  is  measured.  The  Brinell  hardness  is  then  given  by  the 
ratio  of  the  force  P  to  the  curved  contact  area  of  the  defor¬ 
mation  : 

p 

HB  =  f  (B  -  /D2  :  a5)  (1) 

where  D  is  the  diameter  of  the  sphere  and  d  the  diameter  of 
the  permanent  deformation. 


In  most  cases  the  Brinell  hardness  is  not  a  constant 
but  a  function  of  the  load  P  and  the  diameter  of  the  sphere. 
According  to  general  physical  principles  one  would  expect 
that  for  geometrically  similar  deformations,  the  Brinell  hard 
ness  number  is  constant  and  independent  of  their  absolute  mag 
nitude.  This  was  indeed  found  to  be  so. 

If,  for  example,  a  sphere  of  diameter  causes  a  defor¬ 
mation  d^  and  a  sphere  of  diameter  D2  causes  a  deformation  d2 
then  the  Brinell  hardness  number  will  be  the  same  provided 
that  the  deformations  are  geometrically  similar,  i.e.  the 
central  angle  <ji  (Fig.  3)  is  the  same.  This  is  the  case  if 


dl  d2 


According  to  Tabor  (1950)  the  Brinell  hardness  is  not  a 
very  satisfactory  quantity  because  the  ratio  of  the  load  P  to 
the  contact  area  does  not  give  the  average  pressure  over  the 
contact  surface. 

If,  on  the  other  hand,  we  use  the  vertical  projection  of 
the  contact  pressure  p  only  and  apply  the  equilibrium  condi¬ 
tion  we  obtain: 

f3  2 

P  =  2 it  i  p  x  d  x  =  irpa 


where  a  is  the  radius  of  the  circular  indentation  at  the  sur¬ 


face  and  the  projected  pressure  p  is  given  by  p 
measure  of  hardness,  namely: 


This 


(2) 
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was  first  proposed  by  Meyer  and  is  generally  known  as  Meyer 
hardness.  Meyer  also  found  an  empirical  relation  between  the 
load  P  and  the  diameter  of  the  indentation  as: 

P  =  kdn 

where  k  and  n  are  constants  of  the  material.  The  exponent  n 
usually  being  larger  than  2.0  lying  in  the  range  2.0  <  n  <  2.5. 
If  the  load  is  applied  by  spheres  of  different  diameters  then 
the  values  k  and  n  change  also  and  we  get 

n  X  fig 

P  =  k1dj  =  k2d2  =  k3d3  =  - - -  (3) 

For  a  large  number  of  tests  Meyer  found  that  the  exponent  n 
is  independent  of  D  but  that  k  decreases  with  increasing 
sphere  diameters  such  that 


A  =  k1D1n"'2  =  k2D2n_2  =  k3D3n"2 

where  A  remains  constant.  Therefore  we 

p  -  A  jn  _  A  n  _  A 
r  "  D  n-2  Qi  ~  D  n-2  a2  ~  D  n-2 
12  3 

from  which  follows  that: 


may 


,n 


write : 


(4) 


(5) 


P_ 

d2 


and 


ri  n-2 

A  (4, 


A  (4)n 


( 6a) 


( 6b) 


Equation  (6a)  states  that  for  geometrically  similar  deforma¬ 
tions  ($  constant)  the  hardness  number  must  be  constant 


10 


whereas  (6b)  indicates  that  for  geometrically  similar  defor- 

P 

mations  the  ratio  — -r  must  be  constant. 

d2 

That  is  to  say  a  static  load  of  3000  kg  applied  by  a 
sphere  of  10  mm  causes  a  plastic  indentation  which  is  geo¬ 
metrically  similar  to  that  caused  by  a  load  of  750  kg  applied 
through  a  sphere  with  D  =  5  mm  or  by  a  load  of  30  kg  with  D  = 
1  mm.  In  all  these  cases  the  hardness  number  must  remain  the 


same  . 


For  the  general  case  one  may  write: 


JL  =  F(i) 

,  2 


(6c) 


where  F  is  a  suitable  function  expressing  the  similarity 
principle . 

In  order  to  avoid  the  problem  of  the  different  geometry 
for  different  penetration  ratios  ^  which  presents  itself  when 
measuring  the  Brinell  or  Meyer  hardness,  P.  Ludwik  (1903) 
chose  a  conical  indenter  with  an  angle  of  90°  instead  of  a 


sphere.  Thereby  the  penetration  t  is  equal  to  the  indenta- 

d  ird2 

tion  radius  The  contact  area  is  then  given  by  A  =  - - 

2  4sin45° 

and  the  Ludwik  hardness  number  can  be  computed  from: 


HL  " 


4  P  sin  45°  _  _  P 

- r -  =  0.898  ■ 

it  d^  d^ 


For  a  homogeneous  material  the  hardness  number  of  Ludwik 
is  thus  independent  of  the  penetration  ratio  and  geometric 
similarity  obtains  for  any  load  P  or  any  penetration  depth  t. 
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Ludwik  used  this  hardness  number  mainly  to  investigate 
the  internal  friction  of  metals  (Ludwik,  1916). 

As  a  consequence  of  these  investigations  by  Brinell, 

Meyer  and  Ludwik,  a  large  number  of  methods  for  measuring 
hardness  were  developed  which  in  detail  vary  mainly  in  the 
absolute  magnitude  and  shape  of  the  indenters  and  in  the  form 
in  which  the  loads  are  applied  (statically  or  dynamically). 
Here  shall  be  mentioned  only  the  methods  of  Vickers  and 
Rockwell.  For  dynamic  load  application  a  series 'of  methods 

were  developed,  e.g.  the  Shore  scleroscope  an^  the  pendulum 

/  / 

by  Herbert  and  that  by  Kuznetsov.  A  detailed  description  of 

/ 

/ 

all  these  methods  is  given  by  Williams  (1942). 

We  are  particularly  interested  here  in  to  what  degree 
the  principles  of  hardness  measurement  have  been  applied  in 
the  field  of  soil  mechanics. 

Hardness  Measurement  in  Soil  Mechanics 

A  surprising  number  of  soil  mechanics  tests  can  be  traced 
back  either  directly  or  at  least  indirectly  to  the  principles 
of  hardness  measurements  just  discussed.  For  example,  as 
direct  descendants  we  may  consider:  the  steel  ball  test  by 
Tsytovich,  the  Swedish  cone  test,  the  Proctor  plasticity 
needle,  the  California  bearing  ratio  test,  the  ram  penetrom¬ 
eter  for  measuring  the  hardness  of  a  snow  surface  developed 
by  CRREL  (Cold  Regions  Research  6  Engineering  Laboratory)  of 


12 


the  U.S.  Army,  the  Waterways  Experiment  Station  cone  pene¬ 
trometer  of  the  U.S.  Army  Corps  of  Engineers  and  the  aerial 
penetrometer  of  the  U.S.  Air  Force.  As  soil  mechanics  tests 
which  have  their  origin  at  least  indirectly  in  the  principles 
of  hardness  measurement  may  be  considered  the  U.S.  standard 
penetrometer  and  the  penetrometers  used  in  the  various  coun¬ 
tries,  e.g.  in  Denmark,  Holland,  Sweden,  Switzerland  and 
Germany.  The  various  penetrometer  points  are  shown  for  exam¬ 
ple  by  Terzaghi  and  Peck  (1948)  and  the  experiences  made  with 
these  various  penetrometers  are  described  in  the  literature. 
These  latter  methods  have  the  purpose  of  measuring  the  shear 
or  penetration  resistance  of  a  soil  at  some  depth  either  in  a 
bore  hole  or  after  the  penetrometer  has  been  pushed  into  the 
ground  already  some  depth. 

These  methods  are,  therefore,  of  little  interest  here. 

We  are  mainly  interested  in  those  methods  and  tests 
which  we  classified  as  direct  descendants  of  the  hardness 
measurement  principles.  In  the  following  we  shall  give  a 
short  description  and  discussion  of  the  various  methods. 

1 .  The  Steel  Bail  Test  by  Tsytovich 

To  investigate  the  consistency  of  cohesive  soils, 
Tsytovich  modified  the  Brinell  test  in  the  following  manner. 

A  steel  ball  of  diameter  D  =  9.5  mm  is  forced  by  a  load  P 
into  the  smoothened  surface  of  the  soil.  The  load  is  chosen 
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such  that  the  penetration  :!t"  is  between  0.5  and  1.0  mm  where 
t  is  measured  with  an  accuracy  cf  0.01  mm.  The  apparatus  is 
shown  in  Fig.  4.  The  penetration  is  read  after  a  load  period 
of  5,  10,  30,  60  and  300  seconds  and  the  corresponding  shear 
resistance  is  computed  by  the  formula: 

S  s  0.057  i-  (8) 

Dt 

where  P  is  the  load  in  kg,  D  the  diameter  of  the  sphere  (cm) 
and  t  the  penetration  (cm).  This  formula  was  derived  for 
ideal,  plastic,  i.e.,  frictionless  materials,  and  is  supposed 
to  be  usable  for  soils  with  a  friction  angle  of  up  to  7°. 

For  soils  with  a  higher  friction  angle  the  results  must  be 
multiplied  by  a  correction  factor  m  according  to  Beresanzev* 

friction  angle  correction  factor  m 

10°  0.61 

20°  0.28 

30°  0.11 

Since  the  penetration  depth  t  is  time-dependent  and 
increases  with  increasing  load  duration  this  method  yields 
a  shear  resistance  that  decreases  logarithmically  with  the 
load  duration. 

This,  indeed,  also  corresponds  to  actual  clay  soil 
behavior  which  exhibits  a  shear  strength  that  tends  assymp- 
totically  towards  an  ultimate  value  as  the  duration  of  the 
load  is  increased.  It  is  interesting  to  note  that  this 
method  is  one  of  the  few  ones  which  is  able  to  reveal  the 
rheological  properties  of  a  soil. 
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Nevertheless,  a  certain  amount  of  the  penetration  t  may 
also  be  caused  by  consolidation  which  will  introduce  a  small 
error;  however,  its  magnitude  can  be  determined  fairly 'well, 
if  it  should  be  considered  necessary. 

2 .  The  Swedish  Cone  Test 

This  test  follows  exactly  the  same  principle  as  that 
of  Tsytcvich  except  that  the  sphere  is  replaced  by  a  cone 
with  an  angle  of  60°  or  30°  respectively.  The  apparatus  is 
shown  in  Fig.  5.  The  disturbed  or  undisturbed  soil  is  placed 
in  the  container  and  the  surface  is  leveled  and  smoothed. 

The  cone  is  then  adjusted  by  a  screw  such  that  its  point  just 
touches  the  soil  surface.  Then  the  cone  is  released  and  pene¬ 
trates  by  its  own  weight  into  the  soil.  This  penetration  is 
measured  by  a  micrometer  and  serves  to  determine  the  strength 
index  H  of  the  soil  which  is  defined  as  one  sixth  of  the  load 
producing  a  penetration  of  1  cm  of  the  60°  cone. 

From  a  large  number  of  tests  a  Table  could  be  established 
giving  the  strength  index  for  any  penetration  and  load.  The 
strength  index  H  is  a  purely  empirical  value  and  serves  mainly 
to  determine  the  sensitivity  of  a  clay  and  to  compare  the 
relative  strength  ci  various  soil 

If  and  H2  is  the  strength 
and  disturbed  state  respectively. 


layers . 

index  in  the  undisturbed 
then  the  sensitivity  is 
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H.,  /Ho 

.  In  the  range  of  small 

in die i  an  empirical 

relat ion  - 

ship 

was  found  between  H  and  the  shear  strength  s , 

namely : 

H  tX.  /  2v 

s  a  4oq-  (kg/cm  ) 

This 

relationship  is  supposed 

to  be  valid 

for  H  < 

240  . 

Although  the  experiment  is  evaluated 

on  a  purely  empir- 

ical 

basis  it  is  claimed  to  be 

very  useful 

• 

Table 

1 

Strength  Index  for 

Swedish  Cone 

Test 

Penetration  (mm) 

Strength  Index  H 

load  P 

100  grams 

60  grams 

10  grams 

cone  angle 

30° 

60° 

60° 

2.0 

1140 

195 

- 

3.0 

526 

100 

- 

4.0 

296 

58 

- 

5.0 

189 

36.5 

- 

6.0 

140 

26.8 

4.5 

7  .0 

104 

19.7 

3.3 

8.0 

79 

15.4 

2.6 

9.0 

62 

12.2 

2.0 

10.0 

42 

10.0 

1.7 

11.0 

8.2 

1.4 

12.0 

6.9 

1.2 

14.0 

4 ,9 

0 . 8 

16.0 

3.6 

0.6 

18.0 

- 

0.4 

20.0 

- 

CO 

• 

o 

3 .  The  Proctor  Plasticity  Needle 

R.  H.  Proctor  (1933)  proposed  to  measure  the  force 
necessary  to  push  a  plane  indenter  of  circular  cross-section 
with  a  penetration  rate  of  about  0.5  inch  per  second  into  the 
surface  of  a  compacted  soil.  The  apparatus  used  is  shown  in 
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Fig.  6.  Exchangeable  indenters  with  cross-sectional  areas 
varying  from  1.0  to  0.01  square  inch  are  pressed  by  a  rod  and 
handle  assembly  through  a  spring  into  the  soil  and  the  neces¬ 
sary  force  is  read  off.  The  necessary  force  should  be  around 
300  psi.  at  the  optimum  water  content.  Proctor  proposed  this 
test  as  a  construction  control  test  co  determine  whether  the 
required  compaction  had  been  obtained  at  the  permissible 
water  content.  He  considered  a  penetration  resistance  from 
200  to  300  psi.  as  acceptable.  The  test,  however,  has  been 
replaced  more  and  more  by  "in  situ"  density  measurements  and 
is  hardly  used  any  longer. 

4 .  The  California  Bearing  Ratio  Test  (CBR) 

In  1 S40  the  U.S.  Army  Corps  of  Engineers  was  suddenly 
faced  with  the  task  to  develop  a  method  that  would  permit  the 
quick  and  reliable  design  of  airfield  and  highway  pavements. 
After  months  of  intensive  research  the  conclusion  was  reached, 
that  the  method  developed  by  the  California  State  Highway 
Department  was  the  most  suitable.  It  was  therefore  adopted 
for  the  design  of  asphalt  pavements  and  it  permitted  to  take 
advantage  of  the  experience  already  gained  with  this  test  in 
California.  The  CBR  test  is,  in  principle,  a  penetration- 
shear  test  which  determines  a  bearing  strength  modulus.  This 
bearing  strength  modulus  then  determines,  according  to  a 
purely  empirical  relationship,  for  a  given  wheel  load  the 


required  pavement  thickness.  For  the  design  of  airfield  pave¬ 
ments  the  design  curves  of  the  California  State  Highway  Depart¬ 
ment  were  extended  and  confirmed  by  measurements  on  existing 
airfield  pavements.  This  method  then  was  used  extensively 
during  and  aft^r  the  second  world  war  and  probably  is  one  of 
the  most  successful  methods  of  modern  soil  mechanics  . 

The  equipment  for  the  CBR  test  is  shown  in  Fig.  7.  In 
this  test  a  circular  steel  cylinder  of  1.954-inch  diameter  is 
pressed  into  the  soil  surface  by  a  hydraulic  piston  and  the 
load-penetration  diagram  is  found.  In  the  normal  case  the 
stress  corresponding  to  the  penetration  of  0.1  inch  divided 
by  the  stress  corresponding  to  the  same  penetration  (0.1  inch) 

for  the  standard  material  (crushed  rock)  is  the  CBR  value: 

0soil 

CBR  =  - - .  Since  the  stress  for  the  standard  material 

^standard 

at  0.1  inch  penetration  is  1000  psi.,  a  soil  stress  of  300 
psi.  at  0.1  inch  penetration  would  give  a  CBR  value  of  30%. 

This  method  has  found  wide  application  in  the  western 
hemisphere  and  is  being  used  with  great  success.  A  detailed 
investigation  of  the  method  was  carried  out  in  1945  by  the 
Waterways  Experiment  Station  in  which  the  soil  deformation 
chare eristics  below  the  piston  were  carefully  examined  also 
(see  Fig.  8).  The  CBR  test  is  thus  surely  that  test  method 
in  which  the  principle  of  hardness  measurement  has  found  its 
widest  application  in  modern  soil  mechanics. 
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5 ♦  The  Ram  Penetrometer 

This  instrument  is  a  steel  rod  with  a  conical  point 
that  is  driven  into  a  material  surface  by  the  repeated  impacts 
of  a  falling  weight  and  the  penetration  is  then  related  to 
the  impact  energy.  In  principle  this  method  is  therefore 
related  to  the  pile  driving  formulae  .  There  are  a  number  of 
ram  penetrometers  in  existence  and  as  a  typical  example  we 
may  just  mention  that  developed  by  the  U.S.  Army  Cold  Regions 
Research  ard  Engineering  Laboratory. 

The  equipment  is  shown  in  Fig.  9.  A  steel  block  of 
weight  W  falls  through  a  height  H  onto  an  anvil.  The  result¬ 
ing  impact  drives  the  conical  point  of  the  rod  into  the  mater¬ 
ial  surface  and  causes  a  penetration  s.  A  detail  of  the  con¬ 
ical  point  is  shown  in  Fig.  10. 

The  ram  hardness  R  is  then  determined  from: 


ww 

R  =  —  +  W  +  Q  (9a) 

s 

where  Q  is  the  weight  of  the  penetromete-  rod.  Since  it  is 
often  difficult  to  determine  the  pene+  nation  s  for  a  single 
blow,  one  often  determines  the  tot-„l  set  from  a  heat  of  n 
blows  and  finds 


(9b) 
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where  sn  now  is  the  penetration  after  n  blows.  The  results 
obtained  with  the  penetrometer  shown  were  compared  with  uncon¬ 
fined  compression  tests  which  yielded  a  linear  relationship 
between  the  unconfined  compression  test  and  the  logarithm  of 
R  according  to: 

Qu  =  4.078  In  R  -  14.72  (kg/cm2)  (10) 

A  series  of  experiments  and  the  scattering  of  the  results  are 
shown  in  Fig .  11  . 

It  is  interesting  to  note  that  in  this  case  the  load  is 
applied  by  impact,  i.e.  dynamically.  As  one  can  see,  equa¬ 
tions  (9a)  and  (9b)  have  the  same  structure  as  the  various 
pile  driving  formulae  and,  therefore,  are  subject  to  the  same 
criticisms  and  reservations  that  were  made  with  respect  to 
those  by  Cummings.  Especially  the  assumption  that  the  pene¬ 
tration  resistance  remains  constant  throughout  the  penetra¬ 
tion  process  along  the  distance  s  is  of  questionable  validity 
and  the  further  assumption  that  the  various  energy  losses  may 
be  neglected  shows  that  these  formulas  at  best  have  an  approx¬ 
imate  value . 


6 .  The  U.S.  Army  Corps  of  Engineers  Cone  Penetrometer 


This  instrument  was  developed  by  the  Waterways  Exper¬ 
iment  Station  of  the  U.S.  Army  Corps  of  Engineers  in  Vicksburg, 
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Mississippi  and  has  been  used  for  several  years  to  determine 

the  shear  strength  or  the  traf f icability  of  a  soil  surface 

for  given  cross-country  type  vehicles.  The  apparatus  was 

already  shown  in  Fig.  2  and  consists  mainly  of  a  conical 

o  2 

point  having  an  opening  angle  of  30  and  an  area  of  0.5  in 
connected  by  a  steel  rod  and  proving  ring  to  a  handle.  The 
cone  is  forced  by  hand  into  the  soil  surface  and  the  force 
required  is  read  off  from  the  proving  ring  dial.  This  dial 
is  calibrated  in  such  a  fashion  that  one  reads  directly  the 
ultimate  bearing  pressure.  This  figure  is  called  the  "cone 
index."  Hechtl  and  Herbst  (1963)  described  the  application 
of  this  index  in  more  detail. 


7 .  The  Aerial  Cone  Penetrometer 

As  we  have  indicated  earlier,  this  instrument  is  the 
logical  follow  or.  to  the  U.S.  Army  Corps  of  Engineers  pene¬ 
trometers.  It  was  developed  particularly  for  those  situations 
where,  for  practical  or  military  reasons,  the  Corps  of  Engi¬ 
neers  penetrometer  cannot  be  used.  The  instrument  (see  Fig. 

1)  essentially  is  a  long  aluminum  tube  having  a  conical  point 
with  a  spike  at  the  lower  end.  The  cone  is  connected  through 
a  spring-firing  pin  assembly  to  a  cartridge  that  sets  off  a 
smoke  signal,  flare  or  color  vane  when  it  is  fired.  If  the 
impact  of  the  penetrometer  is  hard  enough  the  cartridge  is 
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set  off  and  the  resulting  signal  may  be  observed  visually 
from  the  air.  An  infrared  bulb  may  also  be  lighted.  As  the 
spring  of  the  firing  pin  can  be  varied,  one  can  use  stiffen 
and  stiffer  springs  until  one  is  found  which  no  longer  permits 
the  firing  pin  to  fire  the  cartridge . 

The  instrument  has  a  number  of  limiting  characteristics 
and  disadvantages.  For  one,  the  severity  of  the  impact  will 
depend  on  the  impact  velocity.  This  velocity  therefore  must 
be  measured  independently  or,  at  least,  it  must  be  estimated. 
Secondly,  every  impact  cnly  gives  a  relative  indication, 
namely  whether  the  resistance  encountered  by  the  penetrometer 
was  greater  or  less  than  the  one  required  by  the  chosen  spring 
thus  a  determination  of  the  actual  magnitude  requires  at  least 
several,  if  not  a  large  number,  of  impacts  to  bracket  the  true 
value  with  sufficient  accuracy.  One  of  the  most  severe  dis¬ 
advantages,  however,  in  the  practical  application  is  the  aero¬ 
dynamic  instability  of  the  instrument  and  the  geometric 
requirements  for  its  perfect  operation. 

In  order  to  give  perfect  results  the  aerial  penetrometer 
should  impact  upon  the  horizontal  soil  surface  not  only  in 
perfect  vertical  alignment  but  also  without  any  horizontal 
velocity  component.  In  practice,  however,  the  aerial  pene¬ 
trometer  would  most  likely  have  a  small  inclination  a.  Such 
an  inclination  not  only  reduces  any  resisting  force  trans¬ 
mitted  through  the  pin  by  the  factor  cos  a  but  also  causes 
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the  impact  energy  to  be  transmitted  to  a  small  degree  only 
through  the  conical  point  while  most  of  it  is  converted  into 
a  moment  causing  the  penetrometer  to  impact  broadside.  Even 
in  a  perfect  vertical  alignment  impact  such  a  moment  may 
develop,  namely  if  the  penetrometer  has  a  horizontal  velocity 
component  different  from  zero. 

For  the  normal  operational  case  it  will  be  practically 
impossible  to  eliminate  all  these  sources  of  error  and  achieve 
a  perfectly  vertical  impact  with  zero  horizontal  velocity  com¬ 
ponent  and  an  accurately  known  vertical  component. 

The  Princeton  Impact  Penetrometer 

Because  of  these  difficulties  the  author,  in  the  fall  of 
1962,  began  a  research  project  sponsored  by  the  Terrestrial 
Sciences  Laboratory  of  AFCRL  to  investigate  the  feasibility 
of  methods  measuring  the  soil  properties  of  interest  without 
the  disadvantages  cited  above  such  that  aircraft  trafficabil- 
ity  on  a  soil  surface  in  terms  of  predictable  deformations 
could  be  evaluated. 

The  important  practical  as  well  as  theoretical  consider¬ 
ations  for  developing  such  a  method  are  discussed  below. 

Also,  a  description  of  a  laboratory  instrument  (Fig.  12)  to 
check  the  design  concept  is  given  as  well  as  the  results  vf 
several  series  of  tests  that  were  carried  out  with  the  instru¬ 
ment  . 
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The  instrument  itself  was  designed  by  the  author  and 
manufactured  in  the  central  machine  shop  of  the  School  of 
Engineering  and  Applied  Science  of  Princeton  University.  Mr. 
Christian  Hechtl9  research  assistant,  assembled  and  calibrated 
the  electronic  equipment  and  during  the  summer  and  fall  1963 
carried  out  a  first  pilot  test  series  mostly  on  disturbed, 
artificially  prepared  soils  in  the  laboratory.  These  pilot 
tests  had  the  purpose  to  confirm,  in  principle,  whether  the 
chosen  concept  and  system  was  satisfactory,  to  find  the  most 
favorable  shape  of  .he  nose  for  the  drop  capsule  and  to  find 
the  influence  of  water  content  and  dry  density  on  the  test 
results.  The  results  of  these  pilot  tests  were  presented  in 
February  1954  by  Mr.  Hechtl  as  his  Ph.D.  dissertation. 

The  bulk  of  the  results  presented  in  this  report  were 
carried  out  jointly  by  the  author  and  Dr.  Kechtl  during  the 
summer  of  1964  mostly  on  natural,  undisturbed  soil. 

It  is  the  purpose  of  the  remainder  of  this  report  to 
present  first  the  basic,  theoretical  considerations  for  this 
type  of  remote  testing  and  then  analyze  the  test  results  and 
discuss  the  possible  further  improvement  and  development  of 
the  impact  penetrometer  and  its  operation. 


I I J .  The  Drop  Capsule  Penetrometer 
1 .  Test  Principle 

One  of  the  main  difficulties  encountered  with  the 
aerial  penetrometer  is  measuring,  or  at  least  estimating,  the 
impact  velocity  since  the  square  product  of  this  quantity 
determines  the  impact  energy.  A  small  error  in  determining 
this  velocity,  therefore,  causes  a  rather  large  error  in  the 
kinetic  energy  at  impact. 

An  attempt,  therefore,  was  made  to  consider  techniques 
which  would  permit  a  direct  measurement  of  the  velocity. 

This  would  have  been  possible  perhaps  with  microwave  radar 
techniques  using  the  Doppler  principle.  However,  it  was 
feared  that  this  would  lead  to  rather  complicated  electronic 
equipment.  On  the  other  hand  modern  solid  state  technology 
had  developed  a  series  of  accelerometers  using  piezo-electric 
crystals  (which  had  already  found  wide  application  in  ordi¬ 
nance  and  rocket  development)  which  were  uniquely  suited  for 
our  purposes. 

If  one  attaches  an  accelerometer  (measuring  the  deceler¬ 
ation  as  a  function  of  time)  to  any  instrument  impacting  on  a 
target,  one  can  apply  Newton’s  law  and  write: 

F(t)  =  ma( t )  (11) 


25 


and  the  accelerometer  signal  variation  with  time  corresponds, 
except  for  the  mass  m,  to  the  variation  with  time  of  the 
force  exerted  by  the  target  material  on  the  instrument  which, 
however,  due  to  a  second  law  by  Newton  must  be  equal  to  the 
force  exerted  upon  the  target  material.  Moreover,  since  the 
final  velocity  of  a  missile  after  impact  without  rebound  is 
equal  to  zero  one  can  obtain  from  the  acceleration-time  sig¬ 
nal  a(t)  not  only  the  velocity  v(t)  but  also  the  distance 
traveled  s(t).  Because, according  to  an  elementary  rule  of 
differential  calculus: 


v(t)  -[ 


a(t)  dt  +  c. 


and 


s  ( t )  =  Jv(  t )  dt  +  c2 


(12) 


(13) 


If  we  designate  the  instant  of  initial  contact  as  t  =  0 , 
the  velocity  at  this  instant  with  v(0)  and  the  total  period 
of  the  force  exch.  *ge  with  T,  then  according  to  (12)  we  may 
write : 


a(t)  dt  +  c^  =  0  +  c^ 


(0)  =  f 

J  o 
/*T 

(T)  =  /  a(t)  dt  +  v( 0 )  = 


thus 


v(  0 )  =  -  f 

Jo 


a(t)  dt  =  -  A  (T) 


( lua ) 


(14b) 


(14c) 


That  is  to  say  the  impact  velocity  can  be  obtained  as  the 
area  under  the  acceleration-time  diagram.  Furthermore,  we 
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observe  that  for  t  =  0  the  penetration  s(0)  must  be  zero, 
i.e.  crj  -  0  and  the  final  penetration  s(T)  must,  hence,  be: 


s(T) 


v(t)  dt 


(15) 


Assuming,  for  the  sake  of  argument,  the  acceleration¬ 
time  signal  could  be  represented  by  the  first  half  of  a  sine- 
curve,  then: 

a(t)  =  am  sin  tt  ^  (16) 


which  would  then  yield: 


and : 


(17) 


(18) 


Once  the  deceleration -time  signal  during  impact  is  known 
one  also  knows  the  variation  of  the  force  with  time.  The 
impact  velocity  as  well  as  the  final  penetration  can  then  be 
easily  obtained  by  integration  cf  the  sighal. 

The  choice  of  an  instrument  that  gives  an  electronic  sig¬ 
nal  corresponding  to  the  deceleration  upon  impact  thus  avoids 
all  those  difficulties  which  stem  from  the  complicated  meas¬ 
urement  of  the  impact  velocity  or  its  erroneous  estimate . 

It  was  mainly  for  this  reason  that  the  piezo-electric 
accelerometer  capsule  was  chosen.  Here  we  have  to  mention. 
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however,  one  difficulty  that  presents  itself,  namely,  that  of 
the  velocity  components.  A  piezo-electric  accelerometer,  of 
course,  gives  only  the  acceleration  along  its  piezo-electric 
axis.  If,  therefore ,  the  instrument  impacts  with  a  velocity 
component  normal  to  this  axis,  the  deceleration  in  this  direc¬ 
tion  does  not  give  any  signal.  For  the  prototype  of  a  piezo¬ 
electric  penetrometer  this  problem  can  be  easily  solved  by 
using  a  triple  accelerometer  with  their  respective  axes 
mounted  along  the  axes  of  a  cartesian  coordinate  system.  The 
maximum  component  can  then  be  found  by  superposition: 


2  2  2 

a  +  a  +  a 
x  y  z 


(19) 


Since,  for  our  laboratory  experiments,  accelerometer 
axis  and  impact  velocity  component  in  most  cases  could  easily 
be  made  to  coincide  we  avoided  the  complication  of  using  three 
accelerometers.  This  was  also  in  the  interest  of  economy 
since  the  instrumentation  per  channel  cost  about  $4 ,000-$5 ,000 
However,  for  the  prototype  of  an  aerial  penetrometer  a  triple 
accelerometer  or  possibly  a  directional  stabilization  would  be 
necessary . 


2.  Instrument  Capsule 


The  instrument  was  designed  as  a  hollow,  circular 


chromium  steel  cylinder  of  75  mm  diameter  having  an 
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exchangeable  tip  which  was  either  conical,  hemi-sp'nericai  or 
flat.  The  accelerometer  is  attached  by  a  thread  to  a  small 
anvil  inside  the  tip  in  such  a  fashion  that  the  axis  of  the 
accelerometer  and  that  of  the  cylinder  coincide.  The  total 
weight  was  1.573  kg.  In  Figs.  12  and  13  a  photograph  and  a 
section  of  the  instrument  respectively  are  shown. 

A  shielded  cable  transmits  the  signal  through  an  ampli¬ 
fier  system  to  the  screen  of  the  oscilloscope.  The  hori- 
zontr 1  beam  velocity  of  the  oscilloscope  could  be  varied 
widely.  Thus,  the  time  axis  of  the  signal  could  be  arbi¬ 
trarily  extended  or  compressed  which  permitted  a  pioper 
choice  for  the  best  readability  of  the  signal. 

For  the  test  series  used  so  far  the  signal  was  photo¬ 
graphed  by  a  polaroid  camera  which  permitted  to  make  a  judg¬ 
ment  within  a  few  minutes  whether  a  test  had  been  successful. 
Of  course,  the  signal  could  have  also  been  stored  on  a  mag¬ 
netic  tape  which  for  a  prototype  operational  penetrometer 
appears  more  suitable  and  practicable. 

The  first  experiments  with  this  instrumentation  were 
carried  out  by  Christian  Hechtl  under  supervision  of  the 
author  in  the  summer  of  1963.  They  served  above  all  to  con¬ 
firm  the  applicability  of  the  testing  principle  ,  to  calibrate 
the  instrument,  to  test  the  equipment  and  collect  first  exper¬ 
ience  with  it  in  the  laboratory  and  to  find  the  most  important 
parameters  that  determine  the  target  response.  Above  all,  we 
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wanted  to  find  the  most  suitable  shape  for  the  tip  of  the 
capsule . 

The  results  were  reported  in  detail  by  Hechtl  in  his  dis¬ 
sertation.  They  were  summarized  by  the  author  in  November 
1963  in  Status  Report  No.  4,  Contract  No.  AF(l9)-628-2427, 
submitted  to  the  Terrestrial  Science  Laboratory,  and  the  ten¬ 
tative  conclusions  were: 

1.  The  drop  penetrometer,  in  principle,  is  suitable  for 
the  determination  of  soii  properties. 

2.  The  magnitude  of  the  maximum  deceleration  permits 
the  determination  of  the  strength  '  f  the  soil  in  the  range  of 
the  impact  velocities  used  (up  to  15  m/sec). 

3.  First  tentative  comparisons  between  the  Hertz  and 
the  Meyer  theories  of  impact  showed  that  the  results  were 
closer  to  that  by  Meyer.  This  comparison  was  made  by  Hechtl 
using  rather  severe  simplifications. 

4.  The  most  important  parameter  for  the  laboratory  tests 
appeared  to  be  the  dry  density  of  the  soil,  a  second  one  the 
water  content. 

5.  For  clay  soils  with  a  higher  than  Proctor  optimum 
water  content  smooth  curves  resulted  with  a  monotonously 
increasing  and  then  decaying  signal. 

In  very  dry  soils  and  in  sand  and  gravel  we  observed  cer¬ 
tain  vibrations  of  very  high  frequency  which  correspond  to 
the  eigen  values  of  the  vibration  of  the  capsule  and  are 
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superimposed  on  the  target  response  signal.  In  the  normal 
case  these  superimposed  high  frequency  oscillations  have  no 
influence  on  and  do  not  prevent  the  evaluation  of  the  response 
signal . 

Over  and  above  these  characteristic  frequencies  excited 
by  firm  targets  there  also  appear  certain  softening  and 
rehardening  phases  during  the  impact  period,  especially  in 
sands.  These  should  be  further  investigated. 

Later  experiences  with  the  drop  penetrometer  -  especially 
outdoors  -  led  to  the  conclusion  that  it  will  be  impossible 
to  stabilize  the  penetrometer  against  any  lateral  forces  or 
movements  in  such  a  fashion  that  the  instrument  will  always 
impact  vertically  and  without  any  horizontal  velocity  com¬ 
ponents.  To  minimize  any  error  due  to  these  influences,  we 
chose  for  all  future  tests  the  hemi-spherical  tip. 

Among  the  thirty  test  series  reported  by  Hechtl  there 
were,  however,  only  fvui-  v;arriea  out  with  hemi-spherical  head, 
namely  the  series  No.  12,  15,  19  and  22.  The  rest  had  mostly 
been  carried  out  with  a  conical  head  because  it  gives  a  defor¬ 
mation  that,  as  already  observed  by  Ludwik ,  always  maintains 
geometric  similarity  for  any  penetration  depth.  This  initially 
appeared  advantageous  for  the  theoretical  evaluation  of  the 
results . 

The  above  mentioned  experience  and  the  consequent  conclu¬ 
sions  led  to  a  second  test  series  in  which  the  hemi-spherical 
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head  was  used  exclusively.  Initially,  the  penetrometer  was 
to  be  tested  with  various  target  materials  which,  in  contrast 
to  most  natural  soil  deposits,  had  uniform,  well  defined  prop¬ 
erties  that  were  either  known  or  could  be  easily  determined. 
For  this  purpose  we  chose  two  greases  of  varying  consistency 
and  two  paraffins  of  different  hardness. 

Further  tests  were  then  carried  out  on  natural  soils. 

Before  we  discuss  these  tests,  however,  it  will  be  advan¬ 
tageous  to  consider  theoretically  the  state  of  stress  and 
deformation  that  occurs  in  the  target  material  during  the 
impact  process. 

We  deal  here  with  a  typical  case  of  central  symmetry  for 
which  the  use  of  a  cylindrical  coordinate  system  is  most  suit¬ 
able.  We  will,  starting  from  the  theory  of  elasticity,  con¬ 
sider  the  interplay  of  forces  and  the  deformation  using  an 
ideal  plasticity  theory.  Finally,  we  will  derive  relation¬ 
ships  for  a  viscoelastic  half-space  under  the  influence  of  a 
dynamic  load  applied  through  a  spherical  surface. 
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IV.  Basic  Theoretical  Considerations 


1.  States  of  Stress  and  Deformation  with  Axial  Symmetry 


In  any  state  of  stress  or  strain  being  symmetrical 
to  a  rotational  axis  one  uses,  as  is  well-known,  the  cylindri¬ 
cal  coordinates  r  and  z.  We  follow  here  essentially  a  nota¬ 
tion  and  treatment  given  by  Nadai  and  consider  a  volume  ele¬ 
ment  with  the  sides  rd<j>,  dr,  dz  (Fig.  14)  and  designate  the 
displacement  components  of  the  point  P  in  the  radial  and  ver¬ 
tical  direction  with  u  and  w  respectively.  The  resulting 
normal  strains  can  then  be  expressed  as  follows: 


F  -  3u  .  r  -  u  _3w_ 

er  3r  *  t  "  r  »  ez  dz 


(20) 


and  the  shear  strain  in  any  vertical  plane  is  given  by: 


3u  3  w 
rz  3z  +  3r 


(21) 


The  volume  change  (dialation)  is  then: 


e  =  e. 


et  +  ez 


3u  u  3  w 
3r  +  r  +  3z 


(22) 


If  we  designate  the  stress  components  in  the  radial,  tangen¬ 
tial  and  axial  direction  with  or,  and  oz  respectively  ar d 
the  corresponding  shear  stress  with  tpz  =  t  the  equilibrium 
conditions  for  the  volume  element  of  Fig.  14  may  be  stated 
(ignoring  body  forces): 


~  *«  +  »1  ,  o 
r  3  z 


(23a) 
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°z  3t  t  _ 

3z  3r  r 


(23b) 


For  an  elastic  material  the  stresses  can  be  computed  from  the 
strains  using  Hooke's  stress-strain  relationships: 


o  =  2  G  (|H.  +  _»£■._) 
r  3r  1  -  2v 


(24a) 


*t  =  2  G  +  r^r> 


(  24b) 


o  *  2  G  (I*  f  — y.£.....) 
z  3 z  1  -  2v 


(  24c) 


*  *  G  <!f  *  ff>  *  G 


( 2  4d ) 


where  G  is  the  shear  modulus  and  v  is  Poisson's  ratio.  If  we 
introduce  the  stress -strain  relations  of  equations  (24)  into 
the  equilibrium  conditions  (equations  23),  we  get: 


(1  -  2v)(Au  -  ii_)  +  i£  =  o 
„2  3r 


(25a) 


(1  -  2  v )  A  w  +  ~  0 


(25b) 


.  3  ^  3  32  . 

where  A  =  ~  -  +  rj-  +  — ~  is  the  Laplace  operator 
3  r  *■  1  0  i  d  z* 

ingly  the  volume  change 


Accord- 


l-2v 


1  -2v 


2 ( 1+v ) G  ^  cr  +  °t  +  =  E  ^  ar  +  ct  +  0 


(26) 


must  automatically  sat.isfy  the  Laplace  equation 


iE  =  ill  .  iS_  +  1 h. 

3r2  r3r  8z2 


-"■irWfWr  '  in— rn'* 


(27) 


*  ^  —  ti 
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By  using  the  relationship  (27),  the  equations  (25)  can  be 
:ransformed  into  the  differential  equations  of  the  displace¬ 
ment  components  u  and  w: 

(A  -  -L)(A  u  -  -2L)  =  o  (28) 

r  ^  r  *- 

A  A  w  =  0  (29) 

Let  us  consider  now  the  problem  of  an  elastic  half-space 
z  >.  0  which  has  on  its  surface  z  =  0  a  given  distribution  of 
the  normal  and  shear  stresses  or  of  the  deformation  components 
u  and  w.  We  restrict  our  considerations  to  the  case  of  sym¬ 
metry  with  respect  to  a  vertical  axis.  A  problem  of  this  type 
can  be  solved  by  various  methods  especially  by  using  potential 
theory.  General  solutions  are  given,  e.g.,  by  Riemann-Weber 
(1910).  However,  we  consider  here  only  some  partial  integrals 
of  equations  (28)  and  (29)  which  can  be  used  to  express  arbi¬ 
trary  functions  of  r,  either  as  infinite  series  or  as  definite 
integrals.  These  functions  can  be  made  to  satisfy  the  pre¬ 
scribed  boundary  conditions  on  the  surface  of  the  half-space 
z  =  0  as  it  is  often  done  analogously  with  trigonometric  ser¬ 
ies  or  Fourier  integrals  . 

A  double  pair  for  each  partial  integral  of  (28)  and  (29) 
for  u  and  w  respectively  would  be: 


1  =  e±az  J1(ar), 

u2  =  z  e4az  J1(ar) 

(30) 

1  -  e±&2  JQ( 8r)  , 

u?  =  z  e*®*  J1(8r) 

(31) 

i - P* 
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where  JQ(&r)  and  J^(ar)  are  Bessel  functions  of  the  first 
kind  and  of  order  zero  and  one  respectively  for  a  real  vari¬ 
able,  a  and  $  being  real  constants.  These  functions  J^Cx) 
and  J^(x)  as  is  well-known  are  the  solutions  of  the  differ- 

a 

ential  equation. 


for  n  =  0  and  n  =  1  respectively  which  are  regular  at  the 
origin  x  =  0 . 

They  may  also  be  expressed  as  infinite  series  which  are 
convergent  for  any  value  of  x: 

J  (x)  =  1  -  --x / 2  ^  +  <‘-x/ ? ) - ^ x/.2.),.  +  .  .  .  . 

(1  I)2  (2  t)2  (3  ! ) 2 

j  (x>  .  I  -  -  Ol/IL7  +  .  .  .  . 

1  2  i»2i  2 : 3 1  3 : 4 : 


!>**<*• 
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J  Jq(x)  dx 


J.  (x) 
1 


dx 


Jj.Cx)  dx  =  1  ; 


1 


2  .  Surface  Stresses  on  the  Elastic  Half-Space 

We  consider  the  half-space  z  >  0  laving  stresses 
a  =  p(r)  applied  on  its  surface  z  =  0  which  are  symmetrical 
with  respect  to  the  vertical  axis  v  -  *  with  the  shear 
stresses  at  z  =  0  being  zero. 

This  includes  also  the  case  where  the  normal  pressures 
p(r)  act  within  the  circle  r  <  a  and  vanish  for  r  >  a. 

For  such  a  case  the  stress  components  or ,  at ,  oz  and  t 
must  all  decay  to  zero  for  large  values  of  z.  Therefore,  the 
values  in  the  vicinity  of  the  loaded  circle  wculd  be  of  main 
interest . 

We  introduce  the  displacement  components  cf  the  half¬ 
space  : 

I 

u  ~  (A  +  Baz)  e~az  J^(ar)  (36a)  j 


w  =  (C  +  Daz)  e  J0(ar) 


(3  6b 


where  the  values  A-D  are  constants  of  integration  and  u  and  w 
are  solutions  for  (28)  and  (29)  respectively.  The  constants 
A-D  are,  however,  not  independent  of  each  other  but  must  be 
chosen  in  such  a  way  that  the  conditions  (25)  are  satisfied 
simultaneously .  Besides  these  conditions  the  shear  stresses 
on  the  surface  must  also  vanish  which  means  that  for  z  =  G 


+  2”  =  0 
3z  3r 


(  37) 


If  we  carry  out  the  differentiation  of  the  displacements  (36) 
according  to  (37)  we  obtain 


C  =  -  A  +  B 


(38) 


If  we  further  carry  out  the  operations  with  u  and  w 
which  are  required  by  equations  (25)  we  obtain 


B  =  D  =  - 


1  -  2v  * 


(  39a 


.  _  2(1  -  vj 

C  1  -  2v  A 


(39b 


This  means  that  except  for  the  constant  A  all  other  constants 
are  determined  and  the  displacement  components  (36)  can  be 
given  as: 


u  =  A  (1  -  2v  -  az)  e  J^(ar) 


(40a 


w  =  -A  2(1  -  v)  +  oz  e~az  vJ  ( cir ) 


(40b 
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The  constant  A  must  now  be  found  by  a  second  boundary  condi¬ 
tion,  namely  either  by  prescribing  the  distribution  of  the 
normal  stresses  or  that  of  the  displacements  u  and  w  at  the 
surface.  Let  us  assume  that  A  is  a  function  of  a  parameter 
a:  A  (a).  If  we  multiply  the  right-hand  side  of  equations 

(40)  by  da  and  integrate  from  a  =  0  to  a  =  00  we  can  express 
the  displacements  by  the  definite  integrals: 


u  s  f  A(a 


)  (1  -  2v  -  az.)  e~  J^Car)  da 


(  41a) 


CD 

w  =  J*  A  (a)  (2(1  -  v)  +  az]  e 


J  (ar)  da 
o 


(41b) 


which  yields  the  displacements  at  the  surface  z  =  0: 


u  =  (1  -  2v ) 


[  A( a)  J 
Jo 


(ar)  da 


w  -  -2(1  -  v )  j  A(a)  JQ(ar)  da 


( 42a ) 


(42b) 


By  using  the  expressions  (41)  we  can  again  compute  the 


volume  change : 


=  T?  *  7  +  If  =  2(1-2'')  /”*(“)  e'“Z  J0(«r)cd« 


(43a) 


which  for  z  =  0  yields: 


t  ~  2(1  -  2v)  /  A(a)  J  (ar)ada 


(43b) 
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To  evaluate  A(a)  we  consider  the  case  where  the  normal  stresses 
o,  are  a  prescribed  function  of  r  and  apply  the  boundary  con¬ 
ditions  for  z  ~  0. 

According  to  Hooke's  law: 


3w 
3  z 


ve 


1  -  2  v 


(24c) 


we  get 


6  =  -  P  =  f ( r )  -  2G 

z* 


r  A<a>  jo 

Jo 


( ar ) osdot 


A  comparison  of  (43b)  and  (44)  shows  that: 


1  -  2  v 

e  = - - -  p 


(44) 


(45) 


This  means  that  the  volume  change  is  directly  proportional  to 
the  pressure  p.  According  to  the  theory  of  -  ssel  functions 
an  arbitrary  function  f(r)  may  be  expressed  as  a  Fourier 
integral : 

A  00 

f(r)  =  i  JQ(ar)aa 

Jo 


f 

Jo 


f(X)  J0(aX)XdX 


(46) 


By  comparing  equations  (44)  and  (46)  we  see  that  the  function 
A(a)  is  determined  by  the  definite  integral: 


A(ct) 


p( X)  JQ(aX)XdX 


(47) 


accordingljr  the  displacements  u  and  w  of  the  surface  of  the 
half-space  are  given  by : 
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I 


u  --  -  1 ' ~2~G~-~  [  p( X )  XdX  f  J1(ar)  JQ( 
Jo  Jo  1 


aX  )da 


w  = 


1  -  v 

G  "J 


r 


p  (  X  )  X  d  X 


/: 


JQ(ar)  JQ(aX)da 


(48) 


(49) 


The  interior  integral  of  (48)  is  nothing  but  a  discontinuity 
function  of  the  variable  X  which  may  be  omitted  if  we  replace 
the  upper  limit  of  »  by  r  and  hence  we  get: 


u  =  -  1  T.  2V/’ 

2  G  Jo 


p( X ) XdX 


(50) 


Because 


j: 


J1(ar)  JQ(aX)  =  —  for  o  L  X  L  v 


~  0  for  r  L  X  L 


(51) 


We  introduce  an  average  pressure  p  of  the  pressure  p(r)  within 
the  circle  of  variable  radius  (o  ^  r  ^  <») 


-\C 

r  Jo 


p(X)XdX 


(52) 


The  displacements  of  the  surface  z  =  0  becomes 


1  -  2v  _ 
u  =  "  ~~G —  P  r 


(53) 


The  displacements  w,  however,  cannot  be  determined  so  easily. 
At  the  origin  r  =  z  =  0  we  get,  according  to  (49): 


w. 


p( X )  XdX  f 
Jo 


J  (aX)  da 
o 


(54) 


'  ~r  ■ 


- - 
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<- 


But 


from : 

L 


JQ(aX) 


we  get: 


w. 


dX 


(55) 


Thus,  on  the  surface  the  displacements  u  and  the  volume  change 
are  known.  According  to  (24)  we  may,  therefore,  compute  the 
surface  stresses  and  get: 


°z  ~  ~  p(r)  (56a) 

°r  ~  ~  p(r)  +  j  (1  -  2v )  p  (56b) 


ot  =  -  2  vp(r)  -  |  (1  -  2v)  pr 


(  56c) 


which  are  valid  for  every  point  of  the  surface  .  At  the 
origin  z  =  r  =  0  we  get : 

oz  =  -  p ( o )  ( 56d  ) 

or  =  (Tt  =  -  ~  vl  +  2  v )  p(o)  (  56e ) 


3 .  The  Stress  Distribut  len  after  Hertz 

Hertz  investigated  the  uistribution  of  stresses  and 
displacements  of  two  elastic  spheres  which  are  pressed 
together  by  a  force  P  whose  line  of  action  intersects  the 
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centers  of  both  spheres.  If  the  radius  R2  of  the  second 
sphere  approaches  very  large  values,  we  get  the  problem  of 
a  sphere  pressed  upon  the  half-space .  The  pressure  on  the 
surface  is  then  given  by: 

°z  =  ~  P  =  ‘  Po  "  (r/a)2  (57) 


That  is  to  say  within  the  circle  of  contact  with  radius  a, 
the  pressure  is  given  by  (57)  for  all  values  of  r  >  a  the 
pressure  on  the  surface  vanishes.  If  we  would  plot  these 
pressures  as  ordinates  above  the  corresponding  radii,  we 
would  get  an  ellipsoid  of  revolution.  The  resultant  force 
P  which  is  equal  to  the  volume  of  this  ellipsoid  would  be: 


(58) 


According  to  (52)  the  average  pressure  is 


-j.r 

r?Jo 


p  r  d  r 


(59) 


hence 


■  .d 


F  = 


•Poa 

3r2 


t-'-'S 


3/2 


or  o  £  r  $  a 


(60) 


2p0a‘ 

P  =  - 2" 

3r 


r or  r  >  a 


(61) 


The  vertical  displacements  according  to  Foeppi  (1920)  will  be 


v.  =  I  4  (a  -  K_) 

G  8,2  2a 

cl 


(62) 


The  highest  indentation  occurs  at  the  central  axis  r  =  0: 


The  direct  application  of  these  results  is  not  possible 
because  practically  in  all  experiments  with  the  penetrometer 
the  elastic  limit  is  exceeded.  However,  as  will  oe  shown 
later,  the  Hertz  solutions  can  be  used  when  applying  the 
elastic-viscoelastic  analogy  and  it  is  for  this  reason  that 
the  elastic  solutions  are  discussed  in  such  detail. 


4.  Stresses  and  Deformation  with  Axial  Svmmetry  in  the 


Plastic  Range 


As  already  indicated,  the  main  part  of  the  deforma¬ 
tion  of  the  soil  surface  during  the  penetromete r  test  takes 
place  in  the  plastic  range.  Therefore,  an  investigation  of 
the  deformations  due  to  a  sphere  penetrating  into  a  plastic 
body  is  necessary. 

The  plane  strain  problem  of  a  plastic  deformation  of  a 
half-space  by  a  punch  was  investigated  by  Prandtl  (1921)  and 
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Henckv  (1923).  Hencky  also  considered  an  approximate  solu¬ 
tion  for  the  penetration  problem  with  a  symmetrical  axis  and 
noted  that  the  radii  of  the  slip  lines  varied  only  by  12.5% 
from  those  for  the  problem  of  plane  strain. 

In  his  considerations  Hencky  used  the  assumption  of 
Karman-Kaar,  namely,  that  in  the  plastic  range  two  of  the 
three  principal  stresses  are  equal  and  that  the  third  prin¬ 
cipal  stress  is  different  from  the  two  first  ones  by  the 
moment  2k.  That  is  to  say,  it  differs  by  an  amount  corre- 
sponcing  to  the  yield  stress.  He  thus  used  a  statically 
determined  system  of  stress  as  follows: 

a1  s  2k  +  a2  =  2k  +  a3 

Prandtl  gives  an  ultimate  beari'  •»  pressure  of  ps  = 

2k(l  +  <$  )  whicu  for  the  half  space  with  £  =  n/2  gives  p  = 

5 

2k  x  2. SI.  The  considerations  of  Hencky  yield  a  variable 
pressure  (see  Fig.  15)  with  a  value  of  2k  x  3.33  in  the 
center  and  2k  x  2.52  at  the  edge  of  the  punch.  The  average 
pressure  is  2k  x  2.83. 

Because  the  investigations  by  Prandtl  and  Hencky  were 
made  for  a  rectangular  punch,  one  might  assume  that  for  a 
spherical  penetrometer  they  are  not  applicable.  However,  if 
we  compare  Figures  16a  and  16b  which  show  the  slip  li  •„  pat¬ 
terns  for  a  rectangular  and  for  a  spherical  punch  respec¬ 
tively,  we  recognize  that  this  assumption  is  not  justified. 


It  happens  that  in  both  cases  a  conical  wedge  directly 
below  the  punch  is  forced  into  the  material  and  the  angle  at 
the  apex  of  this  cone  must  be  90°  in  both  cases.  This  wedge 
of  material  is  practically  in  a  hydrostatic  state  of  stress 
and  may  be  considered  as  being  part  of  the  penetrometer  thus 
making  the  actual  body  penetrating  into  the  material  identi¬ 
cal.  The  only  difference  is  that  for  the  cylindrical  punch 
the  loaded  surface  remains  constant  while  for  the  spherical 
indenter  the  loaded  surface  increases  with  increasing  pene¬ 
tration.  Iterson  (1947)  analyzed  the  spherical  indenter  for 
the  hardness  test  in  detail.  Using  the  same  assumptions  as 
Prandtl  he  came  to  the  result  that  the  ultimate  bearing  pres¬ 
sure  for  a  spherical  penetrometer  is  also  ps  =  2k  *  2.57  as 
found  by  Prandtl  for  the  cylindrical  punch.  Iterson  also  car¬ 
ried  out  experiments  to  check  Hencky’s  observation  that  the 
radii  of  the  slip  lines  for  the  three-dimensional  problem 
differ  only  by  12.5%  from  those  of  the  plane  problem  (see 
Fig.  15).  He  found  that  the  plastified  zone  D  at  the  surface 
of  the  material  (Fig.  17)  in  the  plane  strain  problem  was 
equal  to  three  punch  diameters:  D  =  3d.  For  the  three- 
dimensional  problem  he  found,  from  several  series  of  tests 
with  different  loads,  different  sphere  diameters  and  penetra¬ 
tion,  an  average  D/d  =  2.76. 

An  exact,  theoretical  analysis  using  the  hypothesis  of 
Karman-Haar  was  carried  out  by  Ishlinsky  (1944).  He  utilized 
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some  mathematical  results  obtained  earlier  by  Sokolovsky 
which  differs  from  that  shown  in  Fig.  15  only  by  the  fact 
that  the  straight  lines  have  a  very  small  curvature.  More¬ 
over,  Ishlinsky  comes  to  the  same  result,  namely,  that  the 
yield  stress  is  approximately  0.34  to  0.36  times  the  average 
pressure.  This  observation  agrees  also  with  the  statement 
of  Hencky  who  gives  the  yield  stress  as  0.35  Pav  • 

The  solutions  by  Hencky  and  Ishlinsky  have  been  criti¬ 
cized  by  Symonds  (1948)  who  claims  that  the  hypothesis  of  the 
two  principal  stresses  being  equal  have  neither  a  physical 
nor  a  mathematical  justification.  This  may  be  formally  cor¬ 
rect,  however  the  application  of  the  Mises  yield  criterion 
carried  out  by  Symonds  shows  that  real  solutions  exist  only 
in  a  few  special  cases  and  for  the  problem  of  plane  strain 
the  solution  is  identical  with  that  of  Prandt 1-Hencky . 

In  summary  we  may  state  that  according  to  the  theory  of 
plasticity,  a  relationship  can  be  established  between  the 
statical  hardness  test  value  and  the  yield  stress  or  the 
shearing  strength  of  a  material.  Hence,  a  spherical  pene¬ 
trometer  can  be  used,  in  principle,  to  evaluate  the  strength 
of  a  material.  Accordingly,  a  number  of  authors  have  inves¬ 
tigated  the  relationship  between  Brinnel  hardness  and  the 
yield  stress  or  tensile  strength.  For  example,  Foster  (1936) 
gives  the  empirical  relationship  between  tensile  strength  ot 
and  Brinnel  hardness  for  non-ferrous  metals  as: 


14 

hb 

=  -jp  -  1  (in  tons  per  square  inch) 

For  steel  and  steel  alloys  the  same  author  gives  the  yield 
stress  as: 

a  =  0.21  to  0.23  Hc 

y  ° 

The  value  0.22  Hg  is  also  cited  by  Hencky  as  the  one 
found  from  experiments  that  correlates  yield  stress  and 
Brinnel  hardness.  The  discrepancy  from  the  theoretical 
value  (0.35  Hg)  is  explained  by  Hencky  as  follows:  “Accord¬ 
ing  to  the  Brirnel  hardness  test,  the  yield  stress  wou  be 
0.22  times  the  average  bearing  pressure.  According  to  our 
calculations  it  should  be  0.35  times  the  average  pressure. 
However,  we  can  hope  to  find  agreement  with  the  Brinnel  test 
only  if  we  assume  the  penetrated  surface  of  radius  AO  ... 
exactly  as  it  results  after  plastic  equilibrium  obtains.  We, 
therefore,  have  to  solve  first  for  different  penetrations  of 
the  sphere,  the  equations  (14a)  and  (14b)  as  well  as  equa¬ 
tion  (12).  Only  then  a  comparison  with  the  experimental 
results  will  be  possible. 

"If  we  assume  the  penetrating  punch  is  bounded  by  a 
conical  surface,  we  see  from  equation  (15)  that  the  pressure 
below  the  tip  increases  logarithmically  toward  infinity. 
Furthermore,  we  note  that  from  the  sector  0BC  the  part  from 

4>  =  -  45°  to  $  =  0°  is  vanishing  more  and  more  whereby  the 

/dz 

—  changes  its  sign  and  yields  a  positive  contri¬ 
bution  to  the  exterior  (negative)  compressive  stress. 


’’The  relationship  between  yield  stress  and  average  bear¬ 
ing  pressure  (based  on  the  projected  surface),  therefore, 
cannot  be  constant  but  depends  on  the  penetration  and  the 
shape  of  the  penetrating  body.  According  to  equations  (12) 
and  (14)  one  could,  at  least  approximately,  determine  the 
relationship  between  yield  stress  and  average  bearing  pres¬ 
sure  for  any  penetration  and  any  penetrometer  shape  and  would 
then  have  to  compare  this  with  the  experimental  results.  By 
such  a  comparison  one  should  also  be  able  to  determine  the 
influence  of  strain  hardening  which  at  present  cannot  be 
taken  into  account  theoret ically  . . . . " 

How  far  these  statements  suffice  to  explain  the  rather 

considerable  discrepancy  between  theory  and  experiment  we 

may  leave  open  at  this  time.  It  is  obvious,  however,  that 

if  a  significant  strain  hardening  of  the  material,  takes  place 

during  the  plastic  penetration,  the  resulting  hardness  number 

Hg  or,  according  to  Hencky,  the  average  bearing  pressure  pav 

will  be  larger  than  it  would  be  without  strain  hardening. 

This  means  that  the  value  will  become  bigger  if  no  strain 
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hardening  occurs  and  may  well  approach  the  theoretical  value . 

Experiments  with  materials  that  closely  approximate  the 
stress-strain  relationship  of  an  ideal  plastic  body  could 
clarify  this  point. 


■  V  *  Dynamics  of  Load  Application  for  the  Drop-Impact 


Penetrometer 

All  considerations  so  far  were  made  for  static  or  quasi¬ 
static  loads.  However,  the  drop  penetrometer  has  a  basically 
different  load  application  because  the  exchange  of  energy 
between  the  penetrometer  and  the  target  material  takes  place 
during  a  very  short  time  interval .  The  loads  are  therefore 
applied  very  rapidly  and  decay  in  the  same  fashion,  i  .e  .  we 
deal  here  with  a  dynamic,  impact  type  of  load. 

1  *  Stereo-Mechanical  Impact 

The  classical  theory  of  impact  is  principally  based 
on  the  conservation  of  momentum  theorem  for  rigid  body  and 
leads  to  a  simple  mathematical  formulation.  The  application 
of  this  theorem,  however,  does  not  permit  the  determination 
of  any  stresses  of  forces  during  the  impact  process.  One  can 
merely  determine  the  initial  and  final  velocity  components 
and  compare  them  with  the  applied  linear  or  angular  momentum. 
Also,  the  classical  theory  is  unable  to  describe  any  local 
deformations  at  the  point  of  impact  and,  moreover,  assumes 
the  part  of  the  kinetic  energy  that  is  transformed  into  vibra¬ 
tions  to  be  negligible.  This  assumption  is  said  to  be  valid 
for  the  impact  of  two  spheres  (Raleigh,  1906)  or  for  that  of 
a  sphere  with  a  huge  rigid  mass  (Hunter,  1957).  For  an 
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ideally  elastic  impact  a  second  condition  can  be  found  by 
using  the  conservat ion  of  energy  theorem  which  permits  the 
determination  of  the  final  velocities. 

If  the  impact  is  not  completely  elastic  and  permanent 
deformations  do  occur,  one  introduces  a  coefficient  of  resti¬ 
tution  "e.”  This  coefficient  describes  the  degree  of  energy 
exchange  such  that  for  e  =  1.0  the  impact  is  ideal  elastic, 
for  e  =  0  it  is  fully  plastic.  The  value  e  may  also  be  found 
as  the  ratio  of  the  velocity  components  normal  to  the  contact 
surface  before  and  after  impact: 

^f inal 
vinitial 

Frictional  forces  that  may  appear  on  the  contact  surface 

can  be  taken  into  consideration  by  a  dynamic  friction  factor 

/ 

f.  This  factor  is  usually  assumed  purely  on  an  empirical 
basis.  ' 

Let  us  consider  a  mass  point  of  mass  m  moving  with  the 
velocity  V  along  a  path  that  may  be  described  by  the  radius 
r.  The  equation  of  the  motion  can  then  be  written  as: 

••  *  M 

m  r  =  m  v  =  P  (65) 

where  P  is  the  force  causing  a  change  in  velocity.  After 
integrating  once  with  respect  to  time  for  the  time  interval 
T  we  get: 

m(v  -  7  )  =  C  ~ 

Jo 


Pdt  =  Q 


(66) 


X 
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where  the  product  mv  is  the  momentum.  The  time  integral  of 

rT_  _ 

the  force  I  Pdt  =  Q  is  called  the  impulse  A  further  inte- 

Jo 

gration  yields: 


m(r  -  r. 


-  vqT)  .[■ 

Jo 


Qdt  =  Gr 


(67) 


where  G  is  the  mean  value  of  the  integral  Qdt.  Equation  (66) 

is  the  impulse -momentum  theorem.  If  we  determine  the  scalar 
1  _ 

product  of  ^(v  +  VQ)  with  both  sides  of  equation  (66)  we  get: 


1  o  „  Os  _  v  t  vQ 

2  m  (v2  -  vQ2)  =  Q  — 5 - 


(68) 


This  means  that  the  change  in  kinetic  energy  is  equal  to  the 
scalar  product  between  the  impulse  and  the  average  velocity 
before  and  after  the  impact. 

If  we  multiply  equation  (66)  scalarly  with  v  and  note 


that : 


v(v  -  vc)  -  ~  (v2  -  vQ2)  +  ~  (v  -  vQ) 


(69a) 


it  follows  that: 


1  _  o  _  7  1  _  _  2 

T  m  -  v  )  =  y  m  ( v  -  v_) 


(69b) 


For  a  rigid  body  the  last  term  of  equation  (69b)  is  usually 
zero  if  there  is  no  friction  and  no  external  forces,  i  .e . : 


1  2  -2,  1  —  %  2 

2  m  (vo  -  v  )  =  J  ™  (v  *  vQ) 


(70) 


This  is  the  well-known  theorem  of  Carnot. 
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) 


) 


If  we  designate  and  v^  the  velocity  of  the  particle  m^ 
before  and  after  impact  respectively  then  the  principal  equa¬ 
tions  of  the  classical  theory  of  impact  for  perfectly  elastic 
conditions  are: 


Zmivi  =  £  miVi 

and 

Zmi(F^i)  =  Zmi(Fi7i) 


(71a) 


(71b) 


Equation  (71a)  postulates  the  conservation  of  momentum  and 
(71b)  that  of  the  moment  of  momentum.  Hence,  for  the  normal 
impact  of  two  bodies  the  following  equations  apply; 


t  t 

mlvl  +  m2v2  s  mlvi  + 


(71c) 


ml ( V1 ) 2  +  m2  ^  v2  ^ 


2  2 

s  ®lvi  +  r,i  2 v  2 


(7  Id) 


The  terminal  velocities  may  thus  be  determined  from; 

t  2n»2  i  2m^ 

vi  =  vi  -  <n-v2>;  n  =  v2  +  57^7  <vi-v2>  <71 

while  the  magnitude  of  the  force  at  impact  according  to  (66) 
will  be: 


:  t 

P  =  mi(v1"V1)  =  -  m2(v2-v2) 


2m^m2 

(v1-v2> 


2m^m2 

m^+m2 


( 71  f ) 
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2.  Ideally  Elastic  Impact  (Hertz  Theory) 


The  consideration  of  the  contact  stresses  between 
two  elastic  spheres  were  extended  by  Hertz  also  to  the  case 
of  a  dynamic  load.  He  made  the  assumption  that  all  the 
energy  during  impact  is  being  exchanged  twice,  i.e.,  that 
the  energy  which  is  transformed  into  vibrations  of  the  impact¬ 
ing  bodies  may  be  neglected. 

For  the  impact  of  an  elastic  sphere  impinging  vertically 
on  an  elastic  half-space  we  combine  equations  (58)  and  (64) 
and  instead  of  the  radius  of  the  contact  circle  we  introduce 
the  approach  of  the  two  masses  a.  Then  the  load  deformation 
relationship  yields: 


P  =  k0  o' 


(72) 


where  a  is  the  sum  of  the  elastic  deformation  of  the  surface 
of  the  half-space  w^  and  the  shortening  of  the  sphere  radius 
Wq  and  k?  expresses  the  elastic  properties  of  the  sphere  and 
the  half-space.  Thus: 


a  =  Wj  +  Wj 


(73) 


If,  instead  of  the  half-space,  we  first  consider  the 
impact  of  two  spheres  with  the  radii  P.*  and  R2  and  the  masses 
m^  and  m2  respectively  then  the  equation  of  motion  can  be 
written  as  follows: 


P  =  -  m^w^  = 


..  mlm2  .. 

-  m2w2  =  qrj  a 


(74) 


Has.  i  Tikt  .  . 
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From  (66)  follows  that: 


* •  ## 

a  =  +  w2 


r  =  0 
x  =  0 


(75) 


By  combining  equations  (72)  and  (74)  and  integrating  using 
the  initial  conditions  a(0)  =  v(0);  a(0)  =  0  where  v(0)  is 
the  relative  impact  velocity  at  the  beginning  of  contact, 
i.e.,  the  sum  of  the  velocities  of  the  two  spheres:  v(0)  = 
v1(0)  +  v2(0),  one  obtains: 


and 


•• 

a 


m-^+m? 
m^m2  k2 


v(0)2] 


3/2 

a 


kjk2 


klk2 


3/2 


5/2 

a 


(76) 


(77) 


The  maximum  approach  occurs  when  a=  0.  consequently: 


The  approach  a  at  any  time  t  is  given  by: 


t 


_ da _ _ 

v(0)2  -  •—  k^k2  a5/2 


(78) 


(79) 


Since  the  deformation  process  is  assumed  to  be  perfectly 
reversible  the  maximum  approach  occurs  exactly  at  the  mid¬ 
point  of  the  contact  period.  This  contact  period  may  be 
found  by  integrating  equation  (77).  For  this  purpose  we 

introduce  a  dimensionless  variable  Z  =  ~~  and  obtain: 

am 


A 


& 


iivu?  -_s  i* 


.  /*1/2t  .  „  r™  _  da  _ 

“21  dt  -  2  /  „  —  -  —  1  ■  .'. 

Jo  J  o  /v(O)^-'—  k^^cr 


2  am  Z’1  d  Z 

V(0)jo  /l Z 


(80) 


The  solution  of  the  integral  yields: 


. .  *  -liL 

5  r(il>  v“° 


=  2.943 


v(0) 


.943  f  5  I2 
.(o)l/5  [4k1k2j 


(81) 


where  V  is  the  gamma  function. 

Hunter  (1957)  gives  a  simplified  approximation  for  the 
relation  expressed  in  equation  (79)  as: 


a  =  a  sin  l-".58 
” 


(82) 


Thereby  the  force  P  becomes  a  function  varying  with  time 
according  to  a  sine-curve: 

1.140  v (  o )  2  .  1  .068  v (  o )  t 

P  =  - r — sm  -  ;  o  <  t  <  ;; — 7 — r 

klam  am  1 *068  v<°> 


P  =  0  ;  t  > 


1.068  v( o) 


(83) 


If  we  let  the  radius  of  the  second  sphere  become 
infinitely  large,  we  have  the  impact  of  a  sphere  on  a  half¬ 
pace  instead  of  the  impact  of  two  spheres. 
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The  equations  (73)  through  (83)  do  not  change  except 
that  the  impact  parameters  k^,  k2,  a,  a,  un,  and  T  are 
simplified  considerably.  A  tabulation  of  these  parameters 
for  the  two  cases  in  question  yields  the  following  values 
(Goldsmith,  1960). 


Table  2 

Elastic  Parameters  for  the  Hertz  Impact  Theory 


impact  centric  impact  of  two  elastic 
param-  spheres  with  radii  and  R2 
eter 


normal  impact  of  an 
elastic  sphere  with 
radius  R^  on  the  elas- 
tic  half-space _ 


rci  +  m2 

in^niQ 


4  (  1  )fh^2 

3  \6 ^+6 2 / / R^+R/ 


3tr  R1R2 


4  Rl+R2 


(fii+fij)  P 


m- 


^flRl 


JL  -i  fir 

3n  6  j  +  6  2  ^ 


x  Ri  p 


a 


m 


2  Rl+R2  3 it 


R1R2 


4a 


(51+62)P 

.2/5 


§7  = <  wp 


f  2  l2/5r  "»l/5  g.  2 

Jl5frv(o)  (  6^  +  6  2  )mlni2  |^1  +  R2  I  jl5irv(o)  (6 

16(m1+m2)  lRlR2  J  16  /r^ 


1+62)i91 


2/5 


f( 

.53  - 


6i+62 )m^m2 

®1  "*’m2 


2/5 


f  Rl~;  r2  1 

v<  o)R1R2  ) 


1/5 


4.53 


f  **  2/5 

I ( 6^+62 )m^ | 

/R1V  (  oT  ) 


The  values  <5^  and  ©2  designate  the  elastic  constants  for 
the  spheres  1  and  2  respectively  for  the  Hertz  impact  problem 
and  are  determined  by: 


6 


(84) 


whore  E  =  Y01  g's  modulus,  v  r  Poisson's  ratio. 

As  stated  initially  these  relationships  are  valid  if  the 
vibration  energy  is  negligible.  According  to  Hunter  (1957) 
the  ratio  of  vibrational  energy  Ey  to  total  energy  fotf  the 
impact  of  a  sphere  on  the  half-space  according  to  (83)  is 
given  by: 


r  r  3 

Ev  -  2.7  I,  3v(oT 

Eo  gc3l4irflgj 

where  g  is  the  gravitational  acceleration  and 
of  sound  in  the  half-space.  The  vibrational 
is  very  small  as  long  as  the  impac\’  velocity 


(85) 

c  the  velocity 
energy,  hence, 
is  smaller  than 


c  . 


The  second  basic  assumption  for  the  Hertz  impact  problem, 
namely  complete  elasticity  and  reversibility  of  the  deforma¬ 
tion,  is  much  farther  off  the  mark.  If  we  compute  e  .g  .  the 
resulting  impact  force  according  to  (72)  with  the  correspond¬ 
ing  values  of  Table  2,  we  get: 


(86) 


which  according  to  (78)  gives  the  maximum  value  of  the  result¬ 
ing  stresses: 
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Pm 


P 

max 


1/5 

- 2  3  R1_1  [30^(61  +  62)m1v(o>'?3 

3  it  z  ■*• 

(87) 


If  vo  let  a  sphere  drop  on  the  half-space  and  introduce 
e.g.  the  values  of  a  steel  sphere  of  10  mm  diameter  and  com¬ 
pute  the  height  of  drop  which  is  possible  without  the  yield 
stress  being  exceeded  in  the  half-space  material,  we  get  the 
following  values: 


i 


Table  3 

Critical  Drop  Height  for  Elastic  Impact 


half-space  material 

Brinell 
hardness  Hg 

critical  height 
plastic  deformation  (mm) 

soft  copper 

55 

0  .0033 

A-7  carbon  steel 

125 

0.182 

high  strength  steel  alloy 

320 

3 . 5 

It  is  thus  obvious  that  for  most  practical  cases  the  assump¬ 
tion  of  complete  elasticity  and  reversibility  of  the  deforma¬ 
tions  is  far  from  realistic. 


3 .  Impact  on  the  Plastic  Half-Space 

Because  of  the  limited  applicability  of  the  Hertz 
theory  and  equation  (72)  an  attempt  was  made  to  formulate 


relationships  that  would  hold  in  the  plastic  range.  As  in 
the  case  of  the  Hertz  theory  we  begin  again  with  the  static 
load-deformat ion  relationship  and  return  to  the  results 
obtained  by  Meyer  who  found  empirically: 

?  =  k  an  (3a) 

where  a  is  again  the  radius  of  the  permanent  crater  and  k  as 
well  as  n  are  constants  of  the  half-space  material  which, 
however,  as  we  remember,  are  also  dependent  upon  the  radius 
of  the  impacting  sphere.  If  the  crater  is  sufficiently 
shallow,  i.e.,  for  penetration  depths  a  <  R/2,  we  may  replace 
the  crater  radius  by  the  crater  depth  according  to  the  approx¬ 
imation  (see  Fig.  3b): 

a2  =  2Rot  (88) 

Then  equation  (3b)  me}  be  rewritten: 

F  =  ka11  (89) 

where  k  and  n  are  now  different  quantities.  Equation  (89)  is 
valid  until  the  process  of  restitution  starts.  The  latter  is 
assumed  to  be  elastic  and  the  relationship  •  tween  the  force 
and  the  deformation  during  this  restitution  process  is  given 
by 


where  ar  is  the  permanent  crater  depth  and  ara  the  maximum 


i- 
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jL 


60 


penetration.  If  the  restitution  is  assumed  to  proceed  accord¬ 
ing  to  Hertz,  the  exponent  s  assumes  the  value  3/2. 

Equation  (89)  may  be  assumed  exact  if  the  elastic  defor¬ 
mation  is  negligible  and  if,  moreover,  the  plastic  flow  pres¬ 
sure  p0  is  acting  uniformly  over  the  entire  contact  surface 
and  remains  constant  throughout  the  penetration  process  (no 
strain  hardening  occurs).  For  this  case 


k  -  2  n  R  pQ  and  n  =  1.0 


(91) 


This  means  the  load-deformation  relation  corresponds  to 

i 

/ 

that  predicted  by  the  theory  of  plasticity.  In  the  general 
case,  however,  these  idealizations  are  not  satisfied  because 
of  strain  hardening,  the  flow  of  the  displaced  material  dur¬ 
ing  the  impact  process  and,  finally,  because  of  the  particular 
stress -s train  relations  that  obtain  for  the  high  strain  vel¬ 
ocities  during  impact.  For  these  reasons  the  values  k  and  n 
of  equation  (89)  will  in  general  also  depend  on  the  penetra¬ 
tion  velocity  a. 

For  the  plastic  penetration  of  a  rigid  sphere  with  radius 
R  into  a  soft  target  material  as  shown  in  Fig.  3b,  one  can 
establish  the  following  simple  relationship: 


P  =  ir  pa2  =  2  7i  p0  R  a  = 


(92a) 
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The  solution  of  this  simple  differential  equation  with 
the  assumption  that  the  flow  pressure  p0  is  constant  and 
independent  of  u  yields: 


(92c) 


The  maximum  crater  depth  occurs  at  the  end  of  the  impact 
period,  namely  when: 


a  =  ±  v(o)  cos 


Vp 

J  m 


t  =  0 


i.e.  for: 


/  2  TT R P0~' 

J  " 


t  ,  T  — 
2  * 


=  Z.  /  m  ~~T 

2  J  2ttRp0 


Thus  the  maximum  crater  depth  becomes: 


‘max 


=  v(o)  f 


m 


J  2irRp0 

Introducing  (94)  into  (93b)  yields: 


T  =  T 


tr  am 


2  v(o) 


=  1  ,57 


hi 


v(  o) 


(93a) 


(93b) 


(94) 


(95) 


As  a  comparison  the  Hertz  theory  gives  a  contact  period  of: 


T  =  2.943 


m 

ToT 


(81) 


and  a  maximum  approach  of: 

1 2  /  5 


a 


max 


5 v(  o ) 2 1 
4klk2  ) 

» 


(78) 


The  formulations  of  equations  (92)  to  (95)  assume  that 


the  dynamic  flow  pressure  pQ  is  a  constant.  They  are  there¬ 
fore  subject  to  the  same  reservations  that  were  made  already 
earlier . 

4 .  Elastic-Plastic  impact 

Tabor  (1951)  considered  the  impact  of  a  rigid  sphere 
on  an  elastic-plastic  half-space.  Assuming  that  the  sphere 
of  mass  m  falls  through  a  height  h-^  onto  the  half-space  and, 
after  impact,  rebounds  to  the  height  h2  leaving  a  permanent 
crater  with  radius  a,  one  can  calculate  the  total  energy  as 
E-^  =  mgh^  and  the  energy  of  restitution  E2  =  mgh2.  If  one 
assumes,  as  before,  the  vibration  energy  to  be  negligible, 
the  energy  used  to  form  the  crater  E3  must  be: 

E3  =  E1  ~  E2  =  mg  0*1  “  (96) 

But  the  energy  of  restitution  can  be  computed  also  from  the 
theory  of  elasticity  because  it  describes  a  purely  elastic 
process  and  one  may  write: 

3  p2 

E 2  =  mgh2  =  Jo  a  u  ^61  +  62^  (97) 

The  volume  of  the  crater,  in  first  approximation,  is  given  by 


(98) 


63 


Rj  is  here  the  radius  of  curvature  of  the  permanent  crater 
which  because  of  restitution  is  larger  than  the  sphere  radius 
.  If  the  flow  pressure  is  again  assumed  constant  the  energy 
Z3  is  given  by: 


E3  =  Po  v  =  P  =  mg  (h1  -  h2) 


(99) 


Expressing  R2  as  a  function  of  by  using  the  Hertz  relations 
one  gets: 


1.1  i 

R2  '  Ri  "  4  a3 


( 6^  +  S2  ) 


(100) 


which  gives: 


E3  =  Po 


__a_ _ 3_  P_ 

4RX  16  a 


TT  (6.,  +  50) 


(101a) 


The  first  term  is  simply  p0Vs  where  Vs  is  the  apparent  volume 

of  the  crater  which  would  result  if  the  crater  had  the  same 

radius  R^  as  the  sphere  .  The  second  term  is  nothing  else  but 
5 

—  E2  as  a  comparison  with  (97)  reveals  at  once.  Therefore: 


E3  *  PoVa  -  |  e2 


(101b) 


from  which  one  gets  for  the  flow  pressure 


Po  = 


mg  (hrf  h2^ 

— 


(102) 


It  is  obvious  that  for  small  values  of  h2 


mghj 


Po  *  -7 


(103) 
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This  result  was  already  observed  experimentally  by  Marke.l 
(1895)  . 

For  the  derivation  of  equations  (96)  through  (103)  again 
a  constant  flow  pressure  was  assumed.  For  the  case  of  a  vari¬ 
able  flow  pressure  one  can  again  return  to  the  relation 
obtained  by  Meyer: 

p  n-2 

p  =  — r  =  K  d  (6a) 

d2 

where  n  is  a  value  between  2.2  and  2.5.  Then  the  energy  of 
plastic  deformation  would  be: 

e3  =  rrr  pV  (104) 


where  now  p  is  the  average  flow  pressure  at  the  end  of  the 


deformation  process. 


n  +  2 

p  =  - -  mg 


2n  -  1  , 

2 ( n  +  2)  2 

V_ 


(105) 


This  formula  yields  values  about  10  to  12%  larger  than  those 
computed  from  equation  (102). 


The  theory  of  elasticity  and  that  of  plasticity 
determine  the  solution  of  problems  of  stress  or  deformation 
by  assuming  idealized  stress-strain  laws.  However,  most 
engineering  materials  follow  these  laws  only  approximately. 
Moreover,  these  solutions  have  an  additional  deficiency, 
namely  they  completely  ignore  an  additional  dimension  -  time. 

These  theories  assume  that  upon  load  application  a  state 
of  strain  is  attained  instantaneously  and,  thereafter,  with¬ 
out  a  change  of  the  load  no  change  of  stress  or  strain  can 
take  place.  Unfortunately,  real  materials  do  not  behave  that 
way.  Time-dependenx  deformation  (creep)  and,  often  combined 
with  it,  changes  in  stress  (relaxation)  do  take  place  and 
these  occur  in  every  material  under  any  state  of  stress.  In 
some  cases,  fortunately,  stress  relaxation  and  creep  are 
negligible  (as  for  example,  for  steel  at  normal  loads  and 
temperatures).  For  softer  materials  as  e.g.  asphalt,  plastics 
and  especially  for  soils  these  time-dependent  deformations 
however  make  up  the  major  portion  of  the  total  deformation. 

For  this  reason  a  consideration  of  these  time-dependent  pro¬ 
cesses  becomes  necessary. 

The  relatively  young  branch  of  mechanics  that  deals  with 
the  t ime-dependent  deformation  of  solid  bodies  is  rheology. 
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As  the  theory  of  elasticity  rheology  formulates  idealized 
stress-strain-time  laws.  With  these  idealized  laws  of  mater¬ 
ial  behavior  solutions  for  a  given  problem  are  found  which 
now  are  time-dependent. 

To  visualize  the  conditions  prevailing  one  often  works 
with  rheological  models  which  are  built  in  such  a  way  that 
the  equations  describing  their  deformations  have  the  same 
mathematical  structure  as  the  stress-stra.i.n-time  laws  of  the 
material  under  consideration.  As  building  blocks  for  these 
models  the  following  elements  are  used: 


element 

symbol 

deformation  equation 

description 

Hooke 

G 

T  *  .  =  G  >T‘  • 

13  *13 

elastic  body 

Newton 

y 

Tij  *  "  Vij 

viscous  fluid 

St .  Venant 

tij  =  0  ;  <  k 

rigid  plastic 

(Ti-i"k)  -  nf  ..  t..  >  k  body 

"*■  J  *  J  A  J 


Combination  of  the  Hooke  and  Newton  elements  in  parallel 
yields  the  Kelvin  model  with  the  deformation  equation: 


i 3 


*  ij  +  "t  i- 


(106) 


combination  in  series  on  the  other  hand  yields  .he  Maxwell 
model  with  the  deformation  equation: 


T i j  =  l/n  t.. 


+  1/G  T. 


13 


(107) 
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By  combining  Kelvin  or  Maxwell  models  in  parallel  or  in  series 
very  complicated  material  behavior  can  be  described  (see  Fig. 
18).  The  corresponding  deformation  equations  can  be  found  by 
the  proper  combination  of  equations  (106)  and  (107).  In 
general  form  the  deformation  equation  of  these  various  models 
can  be  given  as: 


,  a”  *  dn_1 

n  +  pn-l  ,  n-1 
dt“  dtn  1 


+  + 


Pl  *  Po>  T  = 


(  <ln+l 


,n+l  ,n 

^ _  +  q  __SL  + 

.n+1  Hn 


dt1 


dt 


n 


qt  dt  + 


(108a) 


or 

R(t)  =  S(t) 


(108b) 


where  R  and  S  are  polynomials  of  the  operator  d/dt  with  con¬ 
stant  coefficients.  As  can  be  seen,  derivatives  with  respect 
to  time  appear  in  the  deformation  law  which  means  that  in  the 
solution  of  a  problem  partial  derivatives  with  respect  to  the 
time  as  well  as  with  respect  to  the  coordinates  x,  y,  z  will 
appear . 

For  certain  special  problems,  namely  if  the  time- 
dependence  is  linear,  the  problems  may  be  solved  by  intro¬ 
ducing  a  Laplace  or  Fourier  transformation.  By  this  method 
the  problem  is  "translated”  into  a  simpler  problem  in  the 
Laplace  or  Fourier  domain  which  is  now  time-independent.  The 
solution  is  then  carried  out  in  this  domain  and  is  generally 


,ffe.  ..,..4- 
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much  simpler  but  the  results  later  have  to  be  translated  back 
into  the  real  domain  and  in  this  translation  the  time- 
dependence  reappears.  This  method  is  used  frequently  for  the 
simplification  of  complex  mathematical  problems  (Doetsch, 

1937)  and  was  applied  by  Lee  (1955)  to  the  determination  of 
the  stresses  and  deformations  of  viscoelastic  bodies. 

In  doing  this  Lee  found  that  the  viscoelastic  problem  in 
the  real  domain  always  was  reduced  to  the  purely  elastic  prob¬ 
lem  in  the  Laplace,  domain.  This  means  that  problems  solved  by 
the  theory  of  elasticity  can  be  retranslated  into  viscoelastic 
problems  provided  the  method  of  Lee  is  applicable  in  principle  . 
This  requires  that  the  stress-strain  law  be  linear  and  permit 
superposition  as  follows: 

If  the  stress  (^(t)  yields  the  strain  e^Ct)  and  a2(t)  the 
strain  e2(t)  respectively,  then  the  stress-strain  law  is  said 
to  be  linear  if  and  only  if  upon  application  of  o^(t)  +  o2(t) 
the  strain  at  any  time  will  be  c^Ct)  +  e9(t). 

As  in  the  case  of  the  linearity  postulated  by  Hooke’s 
law  in  the  theory  of  elasticity  real  materials  will  probably 
satisfy  this  linearity  condition  only  approximately.  Never¬ 
theless,  the  elast ic-viscoelastic  analogy  has  been  found  to  be 
a  versatile  and  powerful  tool  in  recent  years  for  the  solution 
of  viscoelastic  problems. 


* 
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2 »  Solutions  by  the  Elastic-Viscoelastic  Analogy 

In  1957  Lee  himself  applied  his  method  to  the  con¬ 
tact  problem  of  Hertz  and  found  for  the  contact  stresses  of 
a  sphere  on  a  half-space  subject  to  a  sudden  force  Q  accord¬ 
ing  tc  (108a) : 


R  [p(r  ,t ) 1  =  S  [/a2  -  r2] 


(109a) 


and 

R  [p ( t ) )  =  jjq-  S ( a 3 ) 


(109b) 


As  one  can  see,  the  operators 
the  deformation  respectively. 
Maxwell  model  equation  (108a) 

R  =  -Do  +  Pi  Tt  >  S 
A  comparison  with  (107)  shows 


R  and  S  act  on  the  load  and  on 
For  a  body  behaving  like  a 
yields  : 

=  Vt  <110) 

1 

that  p0  =  —  and  p.^  =  G. 


For  a  suddenly  applied  load  Q  which  thereafter  remains 
constant,  equations  (109)  yield: 


=  Q  R1Q  (P!  +  p. 


t) 


(111) 


The  solution  follows  after  the  translation  back  into  the 
real  domain  which  is  very  simple  in  this  case: 


p(r,t)  = 


it  G 


'  *!o e'  ®<t’T>-/^j  dT 


(112) 
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The  resulting  sti'ess  distribution  is  compared  with  that  of 
Hertz  in  Fig  .  19  . 

For  the  half- space  behaving  like  a  Kelvin  model  the 
relationship  was  found  in  a  similar  manner  and  gave: 

3  3  Rl«  -  -  t 

a  -  8  —  O  -  e  n  )  (113) 

The  solutions  just  found  still  have  some  important  deficiencies 
which  do  not  permit  their  application  to  the  drop  penetrometer 
problem-  First,  the  Maxwell  or  Kelvin  models  are  still  too 
simple  to  describe  a  somewhat  complex  material  behavior  and, 
secondly,  Lee  assumed  for  his  solution  that  the  load  Q  was 
applied  according  to  a  unit  st<*p  function  and  remained  con¬ 
stant  thereafter.  Both  these  conditions  are  different  in  the 
drop  penetrometer  experiment,  although  one  might  consider  the 
results  (112)  and  (113)  as  rather  elementary  and  rough  approx¬ 
imations  . 

Below  we  will  derive  a  solution  for  a  mors  complete  rheo¬ 
logical  model  for  the  half-space  and  for  a  force  that  increases 
and  decays  monotonously  with  time,  i.e.  one  that  varies  as  a 
sine  curve  according  to  equation  (83). 
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Inspection  of  equations  (106)  and  (107)  shows  at 
once  that  for  a  constant  load  on  a  Maxwell  model,  e.g.  the 
deformations  increase  linearly  with  time  without  any  limit . 
Upon  unloading  there  is  an  immediate  elastic  rebound  which 
corresponds  to  the  instantaneous  deformation  exhibited  when 
the  load  was  first  applied.  The  Kelvin  model  cn  the  other 
hand  does  not  have  any  permanent  deformation  upon  unloading 
and  merely  exhibits  delayed  elasticity.  The  Kelvin  body 
requires,  however,  a  certain  time  until  all  the  deformations 
are  recovered. 

Since  the  mathematical  difficulties  increase  very 
rapidly  with  an  Increasing  number  of  model  elements  -  as  we 
shall  see  shortly  -  it  is  advantageous  to  choose  the  simplest 
model  that  still  is  able  to  represent  the  material  behavior 
adequately  and  with  sufficient  accuracy.  The  simplest  model 
that  satisfactorily  expresses  the  mechanical  behavior  of  a 
soil  is  the  four-parameter  model  (Fig.  20).  It  may  be  con¬ 
sidered  as  a  combination  of  a  Maxwell  and  a  Kelvin  model  in 
series  or  a  model  according  to  Fig.  18c  with  Kelvin  elements 
in  series  and  n  =  3  where  the  first  and  last  element  are 
degenerated . 

The  differential  equation  of  the  four-parameter  model  is: 
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and  its  Laplace  transformation  is  given  by: 


s2x  -  s t ( o )  -  t(o)  +  (  —  +  +  — )  [sT  -  t(o)]  +  —  t 


ni  n  2  n2 


n2 


6-  [s2J  -  sfto)  -  T(o)]  +  G±  (sf  -  f(o)] 

A  x  n  2 


(115) 


If  we  introduce  the  initial  conditions: 

# 

t(o)  =  x(o)  =  0;  hence,  also  fl(o)  =  ^(o)  =  0 
we  get: 


2  A  ,G1  ^  G1  G2 


+  (~±-  +  -i.  +  — )  s  +  —  —It  =  [g-s2  +  s]f  (lib) 

n,  n9  n2  j  [  i  n2  J 


Then  the  viscoelastic  analogy  yields 


T  =  2  p  t 


(117) 


where  the  viscoelastic  analogy  modulus  y  is  given  by: 


V  = 


n  G2 
S  *"  +  — — -  s 
H2 


G  -i 

+  ( JL  + 

nl 


G1  g2 


Go 

4  _£)S  + 

H  2  92  ^1^2 


(118) 


We  now  return 
equation  (62) 

"  =  I  (1 


to  the  elastic  problem  of  Hertz  and  rewrite 
as  follows: 


(119) 


By  introducing  Lame’s  constants  X  and  y: 
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G  =  y 


X 

=  2  (  X  +  y) 


equation  (119)  is  rewritten  as 


_3_  X  t  2p _ P  rf_ 

w  ~  16  y  (  X  +  y)  _ 2  (a  "  2a; 

a 

At  r  =  0 ,  the  total  penetration  of  the  drop  capsule 
obtained  as 


(120) 

can  be 


3  X  t  2n  P 
16  y(X  +  y)  a 


(121) 


For  a  spherical  shaped  head  of  the  drop  capsule 


a  =  /W(R  -  -)' 

4 

It  is  thus  obvious  that  a  non-linear  relationsh 
between  vr  and  P  such  that  the  Laplace  transform 
applied  directly.  However,  we  can  introduce  a 
b  =  wa,  which  varies  linearly  with  the  force  P: 


(122) 

ip  exists 
cannot  be 
quantity,  cay 


b 


wa 


JL  x  1  p 

16  y ( X  +  y)  F 


(123) 


If  P  varies  with  time  as  indicated  by  (c3),  P  =  Pm  sin  ut , 
we  have: 

3  X  +  2y  _  . 

b  =  16  yTT“T~y7  Pm  sin  wt  (124) 

Then  the  viscoelastic  solution  in  the  Laplace  transform 


domain  is 
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“  ,  y  3Pm  X  +  2)1  «. — ; - , 

b(s)  =  06~  5TTTTT  lsln  ut| 


(125) 


Letting 


T  X  +  2  yT 

*  =  pTTTTT 


(126) 


and  introducing 


X  =  K  -  o’  y 


(127) 


we  obtain 


♦  (S)  =  T, 


x  t  2j 7 
jj(I  +  y) 

2  1  ( 


=  cTTTr|<-r,  -jff 


(co^  +  ku^Js  +  u)^o)2 

„  (o..+kci>0  w1w0 

s2  +  (_i__  +  <o0)s  +  -i-£ 


Q 

[((0^+ka)2)s  +  w^(u23  [(l  +  4r)sz  +  (a)^  +  a)2  +  ku)2  +  4rai2)s  +  (i)^0)2] 

0  Ui+kUo  (O-i  (Oo 

s(s  +  co9)[s^  +  (-= - -  +  io9)  s  +  3 


(128) 


where 


^2 

»  ^  =  p  »  W1  =  ~  »  W2  =  n 

Go  A  n  i  ^  4  o 


and  K  is  the  bulk  modulus  introduced  earlier  instead  of  X. 

The  inverse  Laplace  transform  yields: 

<j>(t)  =  q”  T+T"  |  (  1  +  4F  )<5  (  t  )  +  (  1  +  T  )  (  (o^+kcoje  2  +S^e  ^  +82e  )| 


where 


1(1  1  1  (  1  2 
=  2  j  w2+J7f-(o3i'!-k(02)r  ±  2  jtC02+mr(  ui  +  }cw2  >3  ~ 


(129) 


1  +  r  UJ1U2 
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1 

l+r 


1+2F+4T2 

ko)2 

OH 

i+r  + 

(i>2  -u) 3 

1+2T+4T2  , 

ktu2 

(UJ 

i+r  "  + 

U) 2 

104 

a)lW2-Wi(03~k  0)2^3 
tO  jlj.  —  0)3 

a3la,2~a>l'°4  —  ^<t020J4 
003-0)4 


By  using  the  convolution  theorem,  the  inverse  Laplace  trans¬ 
form  of  equation  (125)  gives  the  viscoelastic  solution  in  the 
real  domain: 


b(t)  = 


3? 
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(13  0) 

From  equations  (122)  and  (130),  the  viscoelastic  deformation 
w(t)  can  be  calculated.  Since  the  equation  (330)  involves  a 
complex  combination  of  the  elementary  model  parameters  6^ ,  Gj , 
>  n2>  K,  it  has  to  be  calculated  by  an  electronic  computer. 
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The  principal  difficult}'  is  that  the  material  parameters  do 
not  appear  in  linear  form  or  combinations  only.  Thus  it  is 
net  at  all  certain  that  a  unique  ,  real  solution  exists  and 
it  is  possible  that  several  different  sets  of  parameters 
will  satisfy  the  experimental  data  with  sufficient  accuracy. 


*r: 
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VII*  The  Vibrations  of  the  Half-Space  due  to  Impact 


1 .  Basic  Equations 


For  a;,  isotropic,  homogeneous  and  ideally  elastic 
body  the  equations  of  motion  for  the  propagation  of  an  impact 
wave  can  be  given  by  the  equations  of  motion  (131)  below. 

For  a  soil  the  first  two  conditions  are  satisfied  frequently, 
at  least  in  a  statistical  sense,  whereas  the  latter  (ideal 
elasticity)  surely  is  an  approximation  only.  These  equations 
are  : 


(  A  +  p  ) 


(A  +  p  ) 


(A  +  y) 


9  £ 

_2 

„  o  U 

air  + 

p  7  u  v 

X  =  j 

9t2 

9  e 

_2 

p  92v 

9y  + 

pV  v  + 

Y  =  1 

j 

¥- + 

pV^w  + 

Z  =  < 

P s2? 

9  z 

*  9 12 

(131a) 


(131n) 


(131c) 


where  u,  v>  and  w  are  the  elastic  displacement  components  in 
the  cartesian  coordinate  system  x,  y,  z;  and  e  being  the 
volume  change  as  in  equation  (22): 


9  u  9  v  9  w 
e-9x+9y+3z 


9x2  9y2  9z2 


( 132  ) 


(133) 


is  the  Laplace 
por,?nts  of  the 
Lat.:e  constants 


Operator,  f  the  density 
body  forces  while  A  and 
related  as  follows: 


X 
P , 


Y,  Z  are  the  com¬ 
as  before,  are  the 
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(  Hv  )  ( 1  -  2  v  )  5 


K  -  3  G  ; 


2 ( 1+ v ) 


(134) 


In  an  elastic  material  with  infinite  dimensions  in  all 
directions  a  dynamic  disturbance  causes  two  types  of  waves 
independent  of  each  other  such  that  the  total  displacement  of 
any  point  may  be  considered  as  the  sum  of  the  two  independent 
displacements : 


u  =  u1  +  u2  ;  v  =  v^  +  v2  ;  w  =  +  Wj 


(13  5) 


The  components  of  the  first  wave  are  irrotat ional ,  i.e., 
satisfy  the  conditions: 


3w^  3v-j^  3u-^  3w^  3v^  3u-^  ^ 

3y  3z  3z  3x  3x  3y 


(136) 


which  means  that  they  are  determined  by  potential  functions: 


u  =11-  v  =11-  w  =  H 

U1  3  x  ’1  3y  ’  I  3  z 


(137) 


The  components  of  the  second  wave  are  without  dilation,  i.e.-, 
they  satisfy  the  condition: 


3u2  3  v2  3  v;2 

. .  t  r _  t  ^  _ 


'2  3x  3y  3z 


(138) 


Thus  the  total  displacement s  become: 


u  =  -^1  +  u9  ;  v  =  -|1  +  v9  ;  w  =  +  w2 

3x  2  3  y  2  3z  2 


(139) 


Substituting  these  values  into  the  equations  of  motion  (131) 


yields : 


- 


tJi,a 
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L-t  -  ,2  n2  , 

o  -  a  v  4> 
3 12 


~  2 

3  u2  2  9 

“2"  =  b  7  u2 


where 


3 1 ' 


3  v2  ?  2 

2  =  b  V  V2 
3t  1 

3^w o  2  2 

- -1  =  bZ  V  w9 

3 1 2  2 


2  A  +  2y  2  V. 

a  -  j>  \  b  ~  y 


(  140  ) 


(141a) 


( 141b) 


( 141c  ) 


The  deformations  according  to  equation  (140)  are  pure  volume 
changes  according  to 


i  =  ?  * 


(142) 


For  this  reason  the  wave  given  by  equation  (140)  is  often 
called  dilatational  or  longitudinal  wave.  The  second  type  of 
wave  because  of  equation  (138)  is  free  cf  any  dilatation. 

The  propagation  of  this  wave  yields  only  relative  displace¬ 
ments  whose  components  are  given  by: 


u. 


- 


w, 


3w? 

3  V2 

“  2y 

3  z 

(143a) 

3  u  2 

3  W2 

3z 

3x 

(143b) 

ii2 

3  U2 

(143c) 

3x 

3y 
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For  this  reason  the  wave  given  by  equations  (141)  is 
called  a  shear  wave.  Ir.  a  dilatational  wave  the  movement  of 
a  point  is  in  the  direction  of  wave  propagation,  in  a  shear 
wave  it  is  normal  to  it.  The  equations  (140)  and  (141)  are 
not  coupled  v/hich  means  that  the  two  motions  are  entirely 
independent  from  each  other. 


2 .  Propagation  Velocities 

Integrals  of  equations  (140)  and  (141)  can  be  written 
in  the  following  form: 


—  F  ( t  -  — ) 
r  c 


(144) 


where  r  is  the  radius  vector  of  the  point  under  consideration, 
F  is  an  arbitrary  function  which  satisfies  the  initial  and 
the  boundary  conditions.  The  propagation  velocity  c  =  a  for 
the  longitudinal  waves  and  c  =  b  for  the  shear  waves.  Using 
equation  (134)  these  can  be  compuxed  from: 


-  r  1-v  _  E  _  [  2  ( 1  - v  ) 

' J  (l  +  J)(l-2v)  j>  ~J  (l-2v) 

_  rjuT  _  nr 

~J  2  (  1 +  v  )  f  ~J  J 


(145) 


(146) 


from  which  follows  that  the  propagation  velocities  depend 
only  upon  the  elastic  constants. 


The  ratio  of  the  velocities 

a  _  / 2(1 - v  ) 
b  "J  l-2v 

is,  hence,  always  larger  than  unity.  Barkan  (1948)  gives  the 
following  values  for  these  velocities: 

Table  4 


Wave  Propagation  Velocities  in  Soils 


Soil 

Type 

dens ity 

kg  sec2/cm4  10“^ 

a 

m/sec 

b 

m/sec 

Moist  Clay 

1 . 80 

1500 

150 

Loess  (nat.  water  cont.) 

1  .67 

800 

260 

Dense  Sand  6  Gravel 

1 .70 

480 

250 

Fine  Sand 

1.65 

300 

110 

Medium  Sand 

1.65 

550 

160 

Medium  Gravel 

1 . 80 

750 

180 

3.  Vibrations  Due  to 

a  Point  Load 

Assuming  that  we  deal  with  simple  harmonic  waves 
originating  from  a  point  source  we  will  attempt  to  find  solu¬ 
tions  of  equations  (140)  and  (141)  by  introducing  the  follow¬ 
ing  functions: 

last 

'  e  4>(x,y  ,z)  (147) 


and 
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iojt  iut  iujt 

U2  =  e  U  »  V2  =  S  V’  v;2  =  e  Z 


(148) 


where  w  is  the  exciter  frequency.  Introducing  these  functions 
into  the  wave  equation  the  following  conditions  must  obtain: 


and 


(V2  +  h2)  $  =  0 


(V2  +  k2)  U  =  (V2  +  k2)  V  =  (V2  +  k2)  W  =  0 


( 149 ) 


(150) 


Moreover,  because  of  (138)  the  functions  U,  V,  and  W  must 
also  satisfy: 

*\  T.T 

(151) 


H  +  |V  +  =  o 

3x  3y  3z 


The  solutions  for  equations  (150)  can  be  found  by  the 
substitutions : 


U  - 


3  2 

3x3y  ; 


V  = 


a2  , 

3 

3y  3z  * 


3  2  ij>  2 
W  =  — -■  +  k^ 
3  z2 


(152) 


Introducing  these  values  into  (150)  reveals  that ^  must 
satisfy : 

(V2  +  k2)  \l>  -  0  (153) 


Let  us  now  consider  a  cylindrical  coordinate  system  for 
the  half-space  with  its  origin  at  the  surface  such  that: 

x  =  r  cos  9  ;  y  =  r  sin  9  ;  z  =  z 
for  which  the  Laplace  operator  becomes 


o  2 

2  3  13  3 £ 

=  3r2  +  r  +  3 z2 


(154) 


g?1  yww 


"  V  1 


A  ■% 
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which  yield: 
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F(  O 
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-  2g2  -  k2  P  . 
F(fJ  y  5 


P 

y 


F(C)  =  (2£2  -  k)2  -  4£2  a 8 


Hence,  the  displacement  components  at  the  surface  (z 
be  computed  from  (156): 


Sn 


-C(2C2~k2-2aB) 

F(0 


•J 1  (  C  ,  r ) 


P 

y 


w. 


F(  O 


V$’r) 


P 

y 


(161) 
0)  can 


(162a ) 

(162b) 


Let  us  now  examine  the  case  that  P  have  the  form: 


T.  iwt 

I  e  l .  e  . 

v 

P 

p  =  -  27  5d£ 


Introducing  this  value  into  (162)  and  integrating  we  obtain 
the  displacements  of  the  surface  of  the  half-space: 


v 

2 'fry 


/ 


g2(  2C2-k2-2cx8  ) 
FU) 


JoU’r) 


(163a) 


and 


w 

o 


Pv  /*°°k2 £<* 

2 TTy  /  F(U  JoU’r)  df> 

Jo 


(163b) 


The  integrals  (163)  . annot  be  evaluated  directly  because 
F(0  is  zero  for  certain  values,  e.g.  ,  for  £  -  ±  ^  . 

By  introducing  a  new  variable,  rearranging  and  using 
series  approximation  Shekhter  (1948)  succeeded  in  finding  a 
solution  for  small  values  of  j>=  kr: 


tSBT-t 


*■*  ~*  ^ ' 

.  v 


y 


SfefcSr— *r 
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«o  =  T  eiut  +  i  f2)  (164) 

where  the  functions  and  f2  depend  on  Poisson's  ratio.  For 

v  a  0.5  Shekhter  obtained: 

1  o  5 

f1  =  -  1.00196  j  +  0.0598  J  -  0.00607  J'3  +  0.000243  -  ... 

f2  =  0  .0571  JQ(  1  .047 )  +  0.0474  -  0.00647  j>  2  +  0  .000264  J4 


and  for  v  =  0.25: 

f1  =  -  0.119  |  +  0.0895J-  0.0104  +  0  .00466  j55  -  ... 

f2  =  0.0998  JQ(1. 08777)  +  0.0484  -  0.00595  +  0. 00240 4 

In  general  the  vertical  component  in  the  vicinity  of  the 
exciter  force  can  be  given  in  simplified  form  from  (164): 

w  =  -  AQ  ij)  (f)  sin  (wt  -  t  >  (1C 5a) 

where 


(165b) 


4>(j>)  =  /fx2  +  f22  (  165c) 

and 

X  A 

tan  X  ~  T~  ( 1 6 5 d ) 

2 

The  function  iKj3)  for  Poisson'e  ratios  of  v  =  0,  0.25  and 
0.5  are  shown  in  Fig.  21. 

Sauter  (1950)  evaluated  the  integrals  of  the  displacement 
components  for  the  surface  of  the  half-space  due  to  an 
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4) 

or 


he 


aperiodic  impact  and  obtained  the  following  result  for  a 
purely  normal  impact: 


wD  =  - 


')  K  * 

2*fb  (1  -  ^)t 


(166) 


where  is  the  impact  stress. 

If  we  introduce  the  values  given  by  Barkan  for  dense 
sand  (  j’  =  1.70,  a  =  480  m/sec,  b  =  250  m/sec)  and  for  moist 
clay  (  1.80,  a  =  1500  m/sec,  b  =  150  m/sec),  then  we  get: 

w  =  0.0002  — -  for  sand  (167a) 

o  t 


w„  =  0.0004 


ft 


for  clay 


(167b) 


■  5a ) 

55b) 

5  5c) 

55d) 


tnent 


4 .  Damping,  Dispersion  and  Distortion 

If  there  are  also  frictional  forces  active  in  the 
half-space  in  addition  to  the  elastic  forces  that  tend  to 
return  the  body  to  its  original  shape,  then  damped  vibrations 
result.  If  these  frictional  forces  are  assumed  to  be  propor¬ 
tional  to  the  point  velocity,  the  following  equation  describes 
this  type  of  vibration: 

=  a2  V2^  +  25  (168) 

where  6  is  the  damping  factor  which  is  always  negative  .  For 
the  sake  of  simplicity  we  may  consider  a  one-dimensional 


35^  .TF*"':! 
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vibration  according  to  (168).  It  is  advantageous  to  exchange 
the  variables  by  substituting: 


This  means  that  the  damping  factor  is  a  function  of  the  fre¬ 
quency  from  which  follows  that  different  frequencies  are 
damped  in  a  different  way.  This  phenomenon  is  called  dis¬ 
tort!.  '.  Thus  the  wave  front  is  changing  as  it  progresses. 

Let  us  consider  the  solution: 


u 


(174) 


where  U  is  a  function  of  x  only.  Then, 
( 168  )  ,  we  have : 


according  to  equation 


88 


0  ( 


2  iw6 


(  1 7  5  ) 


This  equation  has  two  solutions  of  which,  however,  only 
one  has  a  practical  significance;  namely,  the  one  that  does 
not  yield  infinite  values  for  u  for  large  values  of  x .  If  we 
put : 


2 

w 


-  2iw6 


(a  +  iB)2 


(176) 


the  solution  becomes: 


u  =  BeiU  *  i6)x 


(177) 


which,  according  to  equation  (174),  yields: 

„  -Bx  i(ax  -  wt) 
u  =  B  e  e 

Since  B  is  also  related  to  the  damping  factor,  both 
as  B  are  functions  of  the  frequency  and  are  relateu 


(178) 

a  as  well 
by  : 


2 

a 


and  aB 


(179) 


When  a  is  a  function  of  the  frequency-  the  phenomenon  is 
called  dispersion,  when  B  is  dependent  on  w,  it  is  called 
distortion.  When  both  phenomena  occur  simultaneously,  the 
wave  energy  is  slowly  dispersed  while  the  shape  of  the 
front  is  continuously  changing. 


wave 


Group  Velocity  and  Phase  Velocity 


where 

these 


where 
pos it i 

u  = 

This  e 
quency 

and 


Considering  two  wave  trains  determined  by: 


32u  _  _ 2  92u 
3t^  " 1  3x^ 


(180a) 


-  2  o  *  2 

3  u  _  2  o  u 

-  a 


3 1 ' 


^  3  x^ 


(180b) 


ai  and  a2  are  two  different  propagation  velocities,  then 
trains  may  be  considered  to  be  of  the  form: 


u1  =  A  cos  [ p ^ ^ —  -  t)  +  E1] 


u9  =  A  cos  [p9(—  -  t)  +  E  9  ] 
z  2  a«  *• 


(181a) 

(181b) 


Pi  P  2 

f1  =  —  >  f2  =  2T  are  two  different  frequencies.  Super- 
on  of  the  two  frequencies  yields: 


u^  +  U2  =  2  A  c 


OS  £  IT 


fl  fl  El+E2 

X  ( -  +  - )  -  H  t  (f1  +  f9)  +  - - - 

a.  2  a2  122 

-  E. 


•  COS 


^TT  X 


(• 


f. 


3  2  a2 


-)  -  TT  t  (f  -  f)  + 


) 


11  -  L2  I 

>  (182) 


quation  describes  two  harmonic  motions,  one  with  a  fre- 
and  propagation  velocity  of: 


f  = 


_  X1 


-  f, 


(183) 


a  = 


al a2  ^  ^1  f  2  ■ 


f la2  t  ^2al 


(184) 


respectively , 


and  the  second  with: 


-  Jfc- 
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and 


a ' 


a^a2 ( f \  -  f 2  ) 
fla2  '  f2al 


(185) 


(186) 


We  note  that  for  the  case  of  =  fj,  equation  (182) 

gives  a  harmonic  motion  that  slowly  changes  amplitudes.  If 

the  maximum  displacement  o £  the  point  x  =  0  occurs  at  time  t , 

then  that  of  the  point  x  occurs  at  time  t  +  and  returns 

r  a  ’ 

with  a  period  —  .  An  interesting  characteristic  of  this 
motion  is  the  continuously  changing  wave  front  and  the  recur¬ 
rence  of  the  maximum  amplitudes  in  certain  intervals. 

If  we  consider  two  motions  with  frequencies  that  lie 
within  a  narrow  band,  say  within  f  and  f  +  Af,  we  have  the 
phenomenon  of  dispersion.  The  \alue  of  a'  of  equation  (186) 
is  the  so-called  group  velocity,  while  the  "a"  given  by  (184) 
is  the  wave-  or  phase- velocity .  To  determine  the  former  we 
introduce  into  (186):  Qj  =  a,  a^  =  a  +  Aa,  f ^  =  f  and 
f^  =  f  +  Af,  and  determine  a*  in  the  limit  as  Af  approaches 
zero.  Then: 


a  *  = 


1  -  f /a 


da 

df 


(187) 


According  to  the  binomial  theorem  and  neglecting  higher  order 
terms  ,  we  get : 


a*  =  a  (1  +  f/a  —  ) 

dr 


(188) 


jj. 


- 


\ r — 


Si 


If  we  introduce  also  the  corresponding  periods,  T,  and 
T  +  AT  and  the  wave  lengths  L  and  L  +  AL,  one  can  easily 
deduce  from  (188)  that: 


a  ’  = 


&<i) 


(189) 


The  group  velocity  is  of  primary  technical  interest 
because  it  determines  the  propagation  of  the  maximum  displace 
ments.  The  situation  is  demonstrated  in  Fig.  23.  There  a 
train  of  waves  is  propagated  like  an  energy  packet  traveling 
with  its  group  velocity  without  essentially  changing  its 
shape . 

However,  because  of  dampening  and  dispersion  the  ampli¬ 
tudes  gradually  become  smaller  and  the  frequencies  lower. 

The  velocities  a  and  a*  thus  become  also  dependent  upon  fre¬ 
quency  or  wave  length.  In  water,  for  example,  long  waves 
travel  faster  than  short  ones  and  the  group  velocity  a'  is 
half  the  phase  velocity  a.  In  elastic  rods,  on  the  other 
hand,  short  waves  run  faster  than  long  ones  and  the  group 
velocity  is  twice  the  phase  velocity.  The  group  velocity 
can  thus  be  larger  as  well  as  smaller  than  the  wave  velocity. 
Moreover,  Lamp  (1904)  has  shown  that  the  two  velocities  a 
and  a’  need  not  run  in  the  same  sense:  the  wave  trains  may 
be  running  toward  the  right  with  such  a  velocity  that  the 
center  of  the  group  moves  toward  the  left. 
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6 .  Impact  with  High  (Ballistic)  Velocities 

If  a  projectile  impacts  on  a  body  with  very  high 
(ballistic)  velocities,  a  crater  develops  whose  volume 
depends  on  the  kinetic  energy  of  the  projectile  and  the 
properties  of  projectile  and  target.  This  is  a  classic 
problem  of  ordinance  research  that  has  been  investigated 
extensively  for  some  time.  The  impact  of  a  meteorite  is  a 
problem  of  a  similar  nature.  Some  recent  results  that 
especially  also  discuss  the  penetration  of  projectiles  into 
sand,  rock  and  sod  are  reported  by  Rinehart  (1954),  Goldsmith 
(1960),  McCarthy  and  Carden  (1962),  and  Kornhauser  (1964). 

Since  our  experiments  were  carried  out  at  relatively 
low  velocities  (maximum  15  m/sec),  the  phenomena  encountered 
at  high  velocity  impact  fall  outside  the  scope  cf  our  con¬ 
siderations  . 
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VIII.  Description  of  Test  Series 

The  idea  of  measuring  the  deceleration  upon  impact  of 
an  instrumented  capsule  by  a  piezo-electric  accelerometer  was 
proposed  by  several  people  independently  at  about  the  same 
time.  Hechtl  (1964),  for  example,  states  in  his  dissertation 
that  such  experiments  were  carried  out  on  sand  and  rock  at 
the  Jet  Propulsion  Laboratory  of  the  California  Institute  of 
Technology.  Their  instrument  was  to  be  used  to  explore  the 
soil  conditions  of  the  lunar  surface. 

However,  Hechtl  based  his  remarks  on  a  preprint  of  a 
paper  by  Dr.  Thormann  which  gave  idealized  response  data.  He 
learned  during  a  personal  visit  in  July  1964  that  actual  tests 
had  not  yet  been  carried  out. 

On  the  other  hand,  McCarthy  and  Carden  (1962)  carried 
out  tests  shooting  an  instrumented  cylinder  through  a  "pneu¬ 
matic  gun"  of  20  feet  in  length  and  6  inches  in  diameter 
against  various  target  materials  such  as  sod,  moss,  wood, 
sand,  lead,  and  concrete  with  velocities  up  to  300  m/sec.  At 
the  same  time  the  author  proposed  independently  to  use  a  simi¬ 
lar  capsule  as  a  new  version  of  an  aerial  penetrometer.  The 
capsule  has  been  described  above  (Chapter  III,  2  and  Figs.  12 
and  13 ) . 

The  first  experiments  were  carried  out  by  Dr.  Hechtl  in 
the  Soil  Mechanics  Laboratory  of  Princeton  University.  Hechtl 
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ran  tnirty  different  series  of  tests  varying  such  experimental 
parameters  as  target  soil,  soil  moisture  content  and  density 
and  shape  of  impact  head. 

The  tests  were  carried  out  by  attaching  the  drop  capsule 
to  an  electrical  release  mechanism  which  was  brougnt  to  the 
desired  height  above  the  soil  target  by  a  pulley-cable  system. 
Upon  impact  the  piezo-electric  signal  was  transmitted  through 
a  trailing,  calibrated  and  isolated  cable,  and  a  cathode  fol¬ 
lower  to  a  Textronic  oscilloscope.  The  oscilloscope  screen 
was  photographed  upon  arrival  of  the  signal  by  a  Polaroid 
camera.  Thus  a  permanent  record  of  the  signal  was  obtained 
that  could  be  evaluated  later.  The  ordinate  of  the  resulting 
photograph  gives  the  deceleration  (to  some  scale  depending 
upon  the  amplification  of  the  scope)  as  a  function  of  the 
time  (shown  on  the  abscissa).  A  schematic  sketch  is  shown  in 
Fig.  22.  The  description  of  the  measuring  components  and 
their  calibration  is  given  by  Hechtl  in  his  dissertation  and 
need  not  be  repeated  here  . 

The  first  experiments  showed  that  a  perfectly  centric 
impact  occurred  rarely  with  a  conical  impact  head  which  led 
us  to  use  a  semispherical  head  for  all  later  tests.  Also, 
the  difficulty  of  keeping  an  oblong  capsule  in  a  vertical 
position  until  impact  under  field  conditions  led  us  to  believe 
that  the  proper  shape  of  a  prototype  capsule  should  be  a 
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sphere.  This  will  avoid  any  problem  of  directional  stabili¬ 
zation. 

The  instrument  was  tested  outdoors  under  field  condi¬ 
tions  impacting  on  natural  undisturbed  soils  and  performed 
well.  The  field  measurements  were  then  carried  out  with  the 
experimental  setup  shown  in  Fig.  22  except  that  the  container 
with  the  soil  target  material  used  in  the  laboratory  is 
replaced  by  the  half-space  in  the  field. 

To  measure  wave  propagation  velocities  of  the  target 
soils  in  the  field  experiments  one,  two  or  (in  test  series 
41)  even  three  accelerometers  were  buried  in  the  soil  at  var¬ 
ious  distances  from  the  impact  point.  They  were  placed  with 
their  sensitive  axis  parallel  to  the  soil  surface  and  directed 
toward  the  impact  point.  Since  the  accelerometers  respond 
only  to  normal  accelerations,  the  resulting  signals  are  from 
waves  of  the  type  described  by  equation  (140).  The  distances 

from  the  impact  point  weie  determined  by  direct  chain  measure- 

\ 

ment  after  impact  and  are  designated  by  r  ,  r2  and  r3. 

Since  according  to  Newton  actio  equals  reactio  w,  the 
penetrometer  signal  gives  not  only  the  force  exerted  by  the 
target  material  on  the  penetrometer  but  also  that  exerted  by 
the  penetrometer  on  the  target.  Thus,  the  signals  of  those 
buried  accelerometers  show  the  decay  and  dispersion  of  the 
stress  waves  resulting  from  an  aperiodic,  dynamic  disturbance. 
The  signals  from  the  buried  accelerometers  were  usually  taken 
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by  a  second  camera  on  a  second  oscilloscope  screen  that  w a s 
triggered  by  the  penetrometer  impact  signa] . 

An  attempt  was  also  made  to  test  the  penetrometer  on 
target  materials  that  had  clearly  defined  and  easily  deter¬ 
mined  mechanical  or  rheological  properties.  This  was  tried 
'  ~ter,  soils  of  different  viscosities,  greases  of  differ- 
■nsistency,  wax  and  paraffin.  Not  all  of  these  experi- 
i..  .s  were  successful.  Dropping  the  penetrometer  into  water  , 
was  tried,  e.g.,  in  the  swimming  pool  in  the  basement  of 
Dillor.  Gymnasium.  Because  of  the  extreme  magnification  of 
the  signal  required  with  the  very  low  impact,  a  basic  60  cps 
frequency  signal  from  the  house  current  was  somehow  induced 
and  superimposed  and  made  the  impact  signal  unrecognizable. 
This  experiment  was  therefore  abandoned.  Oils  as  a  target 
material  also  produced  difficulties.  The  impact  was  often 
so  weak  that  the  electron  beam  of  the  oscilloscope  cathode 
tube  was  not  triggered. 

The  experiments  were  numbered  consecutively.  Continu¬ 
ing  with  the  sequence  begun  by  Hechtl  (30),  thus  they  began 
with  No.  31  and  ended  with  No.  45.  Table  5  lists  the  various 
experiments . 

The  experiments  were  carried  out  jointly  by  the  author 
and  his  assistant.  Dr.  Christian  Hechtl,  during  1963-1964. 
Reproductions  of  the  original  Polaroid  photographs  of  the 
oscilloscope  signals  are  shown  in  the  appendix  (pp.  A-l  to 
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Table  5 

_ List  of  Experiments _ 

Series 

Target  Material  Location  Remarks 


31 

Wax : 

Cambar  M-348  at  31°C 

Princeton  Soil 
Mechanics  Lab  . 

Humble  Oil 
Specs .  Data 
Sheet  DG-36 

32 

Lubricating  grease 
Estan  3  (28°C) 

it  it 

ii  it 

it 

33 

Estan  1 

it  it 

Failed  , 
soft 

too 

34 

Paraffin  (at  28°C) 

it  ii 

35 

Water  (failed) 

36 

Failed 

37 

Silty  sand  at 
soil  surface 

grass  field 
behind 

Princeton  Soil 
Mechanics  Lab . 

natural 

soil 

38 

Silty  sand,  20” 
below  surface 

it  it 

ii 

it 

39 

Uniform  beach  sand 

beach  at  Sea¬ 
side  Heights , 

N .  J . 

natural 

beach 

40 

Failed 

41 

Sandy  silt  and  clay 

McCosh  Circle 
Princeton , N . J . 

natural 

soil 

42 

Sandy  clay 

field  off  U.S. 

1  ,  2  mi  .south 
of  Princeton 

field  plowed 
6  harrowed 

43 

Crushed  rock  3j-l" 

Lambertville , 

N.  J  . 

compacted 

44 

Fine  sand 

Sheppard  * s 

Mill  Dam, 
Bridgeton ,N . J . 

natural 

soil 

45 


Swamp 


'Bridgeton  ,N . J 
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A-32).  The  necessary  calibration  constants  for  evaluating 
these  data  are  also  shown.  The  most  important  results  and 
evaluations  carried  out  so  far  are  offered  in  the  following 
chapter . 

We  have  also  appended  those  four  test  series  carried  out 
by  Hechtl  which  were  conducted  with  a  semispherical  head, 
namely  his  series  No.  12,  15,  19  and  22  (Appendix  B-l  to  B- 
12)  . 

Whenever  possible  standard  reference  tests  were  carried 
out  on  all  target  materials  for  comparison.  For  natural  soil 
surfaces  this  was  the  field  CBR  test.  For  the  laboratory 
soils  we  also  carried  out  CBR  tests  as  well  as  unconfined 
compression  tests.  For  the  wax,  paraffin  and  grease  the 
values  supplied  ny  the  manufacturer’s  data  sheets  were  com¬ 
pared.  Figs.  24  to  27  show  some  of  the  experimental  details. 
For  example.  Fig.  24  shows  the  penetrometer  before  impact  or 
the  paraffin  (Series  34);  Fig.  25  shows  the  penetration  into 
Estan  3  (Series  32);  Fig.  26  shows  the  permanent,  plastic 
deformation  caused  in  Cambar  wax;  and  Fig.  27  shows  the  test 
setup  for  Series  No.  41. 


IX. 


Evaluation  of  Test  Results 


1 .  General  Discussion  of  Test  Data 

The  pclaroid  oscilloscope  photographs  in  many  cases 
(for  example,  in  test  series  12,  22,  31,  37,  38,  41,  42  and 
44)  give  monotonously  rising  and  falling  curves  which  can  be 
approximated  by  a  sine  wave.  The  variatio*  of  the  impact 
force  in  these  cases  can  very  well  be  described  by  a  formula¬ 
tion  according  to  equation  (83).  In  those  cases  where  the 
curves  show  a  significant  variation  from  a  sine  curve  this 
fact  can  be  explained  by  the  particu" *x  test  condition.  In 
test  series  32,  for  example,  the  signal  is  rather  constant 
over  a  considerable  time  interval.  This  was  caused  by  the 
fact  that  in  this  series  the  depth  of  penetration  was  always 
larger  than  D/2.  Thus,  after  an  initial  rise  the  penetration 
resistance  remains  constant  because  the  area  of  penetration 
remains  constant  also  and  the  resistance  encountered  is  that 
of  a  viscous  flow  around  a  cylinder  with  a  semispher ical  head. 
The  last  photo  of  that  series,  moreover,  shows  clearly  the 
penetration  of  the  penetrometer  through  the  entire  pail  of 
grease  (see  Appendix  A-5,  photo  2). 

In  series  34  the  brittleness  of  the  target  material  is 
very  clearly  demonstrated.  While  the  overall  curve  still 
follows  a  sine  curve,  successive  brittle  fractures  cause 


local  deviations . 


T 
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A  significantly  different  signal  from  that  of  a  sine 
curve  is  revealed  for  all  coarse  grained  target  materials 
(test  series  15,  19,  39,  43).  First,  the  electron  beam  is 
broadly  diffused  (compare,  for  example,  Appendix  pages  A-15 
to  A-18,  B-4  to  B-9).  Upon  increasing  the  sweep  speed  of  the 
electron  beam  as  done  in  series  43  (Appendix  pages  A-26  to 
A-28)  one  can  observe  that  this  diffusion  of  the  signal  is 
nothing  but  the  superposition  of  a  high  frequency  vibration 
upon  the  basic  impact  signal  similar  to  the  formulation  of 
equation  (182).  This  high  frequency  wave  signal  can  only 
originate  cither  from  the  natural  frequency  vibrations  of  the 
penetrometer  or  those  of  the  piezo-electric  crystal  induced 
by  a  very  hard  impact.  If  one  considers,  however,  that  these 
vibrations  do  not  affect  the  basic  impact  signal  and  are 
merely  superimposed,  the  correct  impact  signal  may  be  con¬ 
sidered  to  be  the  average  value  within  the  diffuse  band  (see 
page  A-28 ) . 

The  superimposed  natural  frequency  vibrations  can  also 
be  detected  on  less  hard  target  materials  (e.g.,  series  12, 
22,  38,  42  and  44)  although  to  a  much  smaller  extent  and  only 
at  the  beginning. 

Even  more  surprising  than  the  diffusion  is  the  general 
shape  of  the  signal  when  Lhe  target  material  is  sand,  gravel 
or  crushed  stone.  There  exists  a  definite  peak  force  which 
drops  off  rapidly,  after  which  comes  again  an  increase 
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followed  by  a  logarithmic  decay.  This  can  only  be  explained 
in  such  a  way  that  in  these  materials  a  "temporary  lique¬ 
faction"  takes  place  at  the  beginning  of  the  impact  period. 
Since  very  few  grains  are  affected  initially,  the  specific 
energy  per  grain  is  very  large  and  one  can  perhaps  rational¬ 
ize  that  during  this  short  duration  the  grains  partially  lose 
their  frictional  support  and  behave  like  a  quasi-fluid.  Upon 
further  penetration,  however,  the  number  of  soil  grains 
affected  by  the  impact  becomes  quite  large  and  the  available 
energy  is  not  sufficient  to  "liquefy"  all  these  particles. 
Thus  a  stabilization  of  the  penetration  resistance  occurs 
with  a  corresponding  dropping  off  of  the  signal  in  the  final 
phase  of  the  impact  period. 

In  any  event,  from  the  variation  of  the  signals  obtained 
it  is  clear  that  for  the  number  of  tests  performed  so  far  a 
fairly  definite  distinction  can  be  made  whether  the  target 
material  was  a  granular,  "liquef ied-stabilized"  soil  or 
whether  it  was  a  plastic-cohesive  soil,  and  whether  the 
impact  was  hard  or  relatively  soft  depending  on  the  magni¬ 
tude  and  extent  of  the  natural-frequency  vibrations  super¬ 
imposed  on  the  basic  impact  signal. 
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If  the  deceleration  a(t)  is  known,  the  impact  veloc¬ 
ity  can  be  determined  from  equation  (14)  by  integrating  a(t) 
over  the  impact  period. 

Since  the  oscilloscope  photos  -  except  for  some  constant 
-  are  nothing  else  but  a(t)  we  only  have  to  carry  out  the 
integration.  Here  a  small  technical  difficulty  is  encountered 
the  sweep  of  the  electron  beam  must  be  adjusted  in  such  a  way 
that  a  good  resolution  of  the  signal  appears  on  the  oscillo¬ 
scope  screen.  If  the  beam  velocity  is  relatively  small,  a 
high  curve  with  a  narrow  base  results  (e.g.,  photos  of  test 
series  41  with  a  beam  velocity  of  =  0.2  cm/millisecond). 
There  the  signal  can  be  followed  over  a  long  period  of  time. 
Thus  the  end  of  the  impact  period  can  be  clearly  established. 
In  the  limiting  case  of  zero  horizontal  beam  velocity  the 
signal,  of  course,  would  be  only  a  vertical  line.  Thus  it  is 
easily  possible  that  for  a  small  beam  velocity  the  signal  is 
not  sufficiently  resolved  along  the  time  axis  to  permit  inte¬ 
gration  with  sufficient  accuracy. 

If  the  beam  velocity,  on  the  other  hand,  is  relatively 
large  (as  for  example  in  test  series  42  with  =  0.5  cm/ 
millisecond)  ,  the  signal  has  an  excellent  resolution  with 
respect  to  time  but  the  end  of  the  impact  period  dees  no 
longer  appear  on  the  oscilloscope  screen. 
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The  technical  solution  cf  this  dilemma  is  best  solved 
by  recording  the  signal  on  magnetic  tape  rather  than  by  photo¬ 
graphing  an  oscilloscope  screen.  The  tape  can  later  be  run 
at  any  speed  desirable  and  no  part  will  be  lost.  For  an 
operational  penetrometer  magnetic  recording  of  the  signal  is 
probably  essential.  Even  though  such  recording  systems  are 
standard  items,  their  price  is  relatively  high  and  for  cur 
pilot  tests  the  oscilloscope  photos  were  considered  suffi¬ 
cient  . 


The  evaluation  of  the  impact  velocities  from  the  test- 
data  is,  for  this  reason,  not  nearly  as  accurate  as  would  be 
possible  using  che  latest  technology.  Nevertheless,  the 
integration  of  the  oscilloscope  curves  using  a  planimeter 
yields  a  rather  good  agreement  with  the  theoretical  computed 
values.  These  were  computed  by  assuming  the  frictional 


resistance  of  the  penetrometer 

falling  in 
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between  the  theoretically  computed  impact  velocity  and  the  one 
obtained  by  integrating  the  penetrometer  signal.  The  mean 
discrepancy  dm  calculated  from: 
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Table  6  (continued) 


Test 

Series 

Drop 

Height 

v(°)1=ya(t)dt 

v(  °  >  2=/^gh* 

v(o)1/v( o)0 

Discrepancy 

cm 

cm/sec . 

cm/sec . 

% 

44 

275 

556 

522 

1  .061 

+  6.1 

520 

786 

1010 

0 .778 

-22  .2 

520 

713 

1010 

0 .708 

-29.2 

12 

61 

309 

346 

0 .894 

-10  .6 

122 

426 

490 

0 .870 

-13.0 

183 

529 

600 

0 .882 

-11  .8 

305 

655 

773 

0 .845 

-15  .5 

42  5 

713 

911 

0  .782 

-21  .8 

700 

643 

1175 

0.548 

-45  .2 

15 

61 

250 

346 

0 .722 

-27  .8 

122 

331 

490 

0  .677 

-32.3 

183 

368 

609 

0  .614 

-38  .6 

305 

441 

773 

0  .570 

-43  .0 

425 

573 

318 

0.626 

-37 . 4 

19 

61 

233 

346 

0.672 

-32.8 

122 

316 

490 

0  .644 

-35  .6 

183 

279 

600 

0  .465 

-53  .5 

305 

537 

773 

0  .693 

-30.7 

425 

662 

918 

0 .724 

-27.6 

700 

662 

1175 

0.562 

-43  .8 

22 

61 

357 

346 

1.080 

+  8.0 

122 

422 

490 

0 .860 

-14.0 

183 

515 

600 

0.836 

-16 .4 

305 

537 

773 

0  .693 

-30  .7 

425 

654 

911 

0  .718 

-28  .2 
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was  found  to  be  13.8%  for  the  test  series  30  to  45.  This 
discrepancy  has  several  components.  First,  the  direction  of 
impact  and  the  axis  of  the  accelerometers  are  not  always  in 
perfect  alignment.  This  error  can  be  eliminated  as  was  dis¬ 
cussed  earlier  by  using  three  accelerometers  with  mutually 
perpendicular  axes  and  superposition  of  the  signals  according 
to  the  Pythagorean  theorem.  Secondly,  the  drop  height  could 
only  be  maintained  with  an  accuracy  of  ±  2%.  This  was  due  to 
slack  and  elastic  extension  of  the  supporting  cables. 

The  accelerometer  has  a  guaranteed  accuracy  of  0.5%; 
however,  due  to  the  signal  amplification  this  is  magnified 
several  times.  In  addition,  of  course,  the  theoretical  value 
with  which  the  experimental  values  are  compared  are  not  the 
exact  and  true  ones  because  of  air  friction.  Finally,  there 
occurs  a  relatively  large  error  in  the  evaluation  especially 
when  either,  because  of  the  sweep  being  too  fast,  the  end  of 
the  impact  signal  does  not  appear  any  more  on  the  photograph 
or,  because  of  the  sweep  being  too  slow,  the  signal  is  not 
resolved  sufficiently  to  permit  accurate  integration.  These 
possible  errors  can  be  avoided  by  recording  the  signal  on 
magnetic  tape  . 

Increasing  use  and  experience  is  certain  to  lead  to  pro¬ 
gressive  elimination  or  reduction  of  various  sources  of  error 
so  that  the  measurements  will  be  sufficiently  accurate  and 
reliable  for  scientific  purposes.  Here  it  is  perhaps 
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interesting  to  note  that  the  mean  discrepancy  for  the  four 
test  series  by  Hechtl  (No.  12,  15,  19,  22)  computed  in  the 
same  manner  still  have  a  discrepancy  of  31.3%  while  those  of 
our  tests  have  only  13.8%. 

The  author  believes  that  this  is  sufficient  improvement 
so  that  one  may  judge  the  technical  feasibility  of  the  method 
as  positive.  We  believe  that  the  mean  error  can  be  easily 
brought  within  the  range  of  3%  to  5%.  This  is  judged  to  be 
entirely  satisfactory  considering  the  natural  scattering  of 
soil  properties  and  the  general  nature  of  the  problem. 


3 .  Strength  Properties  of  Target  Materials 


As  shown  in  Table  3,  the  critical  drop  height  for 
purely  elastic  impact  is  3.5  mm.  even  for  the  hardest  steels. 
Thus  it  is  superfluous  to  attempt  to  evaluate  elastic  con¬ 
stants  from  the  impact  data. 

The  evaluation  of  the  plastic  flow  pressure  pQ,  however, 
should  give  one  of  the  critical  material  constants. 

According  to  equation  (92a)  the  penetration  during  plas¬ 
tic  impact  in  a  simplified  form  is  given  by: 


where 


a 


v(  o ) 
b 


sin 


bt 


2ttR 


(  92c) 


m 
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The  end  of  the  impact  period  is  determined  by  the  fact 
that  the  penetration  velocity  a  has  decayed  to  zero: 


a  =  v(o)  cos  bT  =  0 


or 


bT  =  -r 


(93b) 


Thus,  the  plastic  flow  pressure  can  be  found: 


IM7 


8R  T‘ 


(191) 


T  signifies  the  duration  of  the  impact.  The  flow  pressure, 
however,  can  also  be  obtained  from  the  maximum  deceleration. 
If  we  differentiate,  for  example,  equation  (92d)  twice  with 
respect  to  time,  we  obtain: 


a  -  -  v(o)b  sin  bt 


(192) 


the  maximum  deceleration  then  is: 


a 


max 


=  -  v(o)b  =  -  v(o) 


f  2ttR  ] 

{—  M 


solving  for  p  we  get: 


m 


2  "R  v  ( 


a  V 

max  j 

v(o)  J 


(19  3) 


(194) 


Thus  the  yield  stress  can  be  evaluated  from  the  test  data 
either  by  equation  (191)  or  by  (194).  If  we  introduce  the 
dimension  of  the  penetrometer  model  as  given  in  Chapter  III,  2 
into  these  equations  we  get: 
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W  1573 
m  =  cr  ~  981 


1.6 


2 

g  sec 
cm 


and 


d  =  0.168 
*  o 


1 

T2 


g/cm‘ 


as  well  as 


Po 


0.068 


a 

max 

v(o) 


g  /  cm" 


(191a) 


(194a) 


It  is  intere  ,  to  note  that  according  to  equation 
(191)  the  yield  stress  is  independent  of  the  drop  height,  the 
impact  velocity  or  the  maximum  deceleration  and  depends  only 
upon  the  impact  duration  T.  The  observed  values  do  indeed 
show  a  surprising  constancy  of  T  for  a  given  target  material. 
Table  7  shows  the  evaluation  of  the  pQ  values  by  both  methods. 
For  a  better  comparison  the  average  values  nQ/'  )r  the 
individual  target  materials  are  compiled  in  Table  8  and  are 
compared  with  the  standard  control  test  results  of  CBR  or  the 
equivalent  unconfined  compression  strength  qu . 

The  values  for  the  flow  pressure  have  an  average  scatter¬ 
ing  of  25.5%  for  pQ  computed  from  equation  (191)  and  of  41% 
for  pQ  computed  from  equation  (194).  This  is  not  too  sur¬ 
prising,  because  equation  (191)  uses  only  one  experimentally 
determined  value  while  (194)  requires  arnax  as  well  as  v(o). 
Thus,  errors  are  being  multiplied.  Furthermore,  not  all  the 
scattering  is  to  be  blamed  on  insufficient  accuracy  of  the 
test  method.  A  certain  scattering  is  to  be  expected  because 
the  very  properties  being  measured  are  not  unique  quantities 
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Table  7 


Determination  of  the  Plastic  flow  Pressure 


Test 

Series 

Drop 

Height 

H 

cm  . 

Impact 

Dura¬ 

tion 

T 

fn  •  s  g  c  * 

1 

2 

T 

10"6secT2 

Po 

( eq .191a) 

kg/cm2 

a 

max 

m/sec2 

v(  o )  | 

1 

cm/sec  . 

[a  1  2 
m 

_  v(  o)^ 

Po 

( eq  .194a) 

kg/cm2 

31 

61 

18 

0 .00309 

0  .52 

570 

346 

26900 

1  .83 

122 

19 

0.00277 

0  .46 

730 

490 

22200 

1  .50 

183 

20 

0.0025 

0 .42 

970 

600 

26200 

1  .78 

30  5 

22.5 

0.00198 

0.33 

1140 

773 

21800 

1  .48 

425 

21.0 

0 .00227 

0  .  37 

1350 

918 

21800 

1  .48 

32 

6 1 

80 

0 .000153 

0  .026 

114 

346 

1090 

0  .074 

122 

85 

0.000138 

0  .023 

148 

490 

900 

0  .061 

183 

90 

0 .000124 

0.021 

160 

600 

713 

0.043 

34 

61 

122 

183 

305 

1  .2 

1  .14 
1.12 
1.12 

0.69 

0 .77 

0  .79 

0.79 

116.0 

129.0 

132.5 

132.5 

5920 

8550 

10250 

13150 

346 

490 

600 

773 

2950000 

3020000 

2920000 

2900000 

201 

205 

198 

197 

37 

275 

8.0 

0.0156 

2.62 

2170 

522 

172000 

11  .7 

275 

8.0 

0.0156 

2.62 

2170 

522 

190000 

12.9 

520 

8  -0 

0  .0156 

2.62 

2850 

1010 

80000 

5.4 

33 

275 

275 

520 

8.8 

9.2 

8 .8 

0.0129 

0.0118 

0.0129 

2.17 

1  .98 

2 .17 

1600 

1  c5  0 

2850 

522 

522 

1010 

94000 

87500 

79500 

6  .4 

5  .9 

5  .4 

39 

275 

45 

0 .000495 

0  .083 

456 

522 

7600 

0.52 

520 

42 

0 .00057 

0  .096 

548 

1010 

29  30 

0  .20 

520 

40 

0.000625 

0.105 

502 

1010 

2480 

0.17 

Ml 

275 

7  .0 

0 .0204 

3.44 

2530 

522 

237000 

15 . 8 

275 

6.5 

0  .0236 

4 .14 

2450 

522 

220000 

15  .0 

275 

7  .0 

0  .0204 

3  .44 

2740 

522 

273000 

18.5 

520 

8  .0 

0.0156 

2.62 

3540 

1010 

122500 

8 . 3 

Table  7  (continued) 


Test 

Drop 

Impact 

Dura¬ 

1 

T2 

Po 

amax  v(o) 

#* 

am 

[2  Po 

Series 

Height 

H 

tion 

T 

10_6sec?2 

(eq.lSla) 

P 

m/sec .  cm/sec  . 

v(  o) 

s.  > 

|  (eq  .194a) 

p 

cm . 

m . se c  . 

kg/cm 

kg/cm 

4  2 

2'  ' 

27. 

52C 

17.2 

3  2.8 

13.2 

0.00338 

0.0061 

0.00575 

0.57 

1.02 

0  .97 

620 

1185 

1510 

522 

522 

1010 

24500 

51000 

2  2  5  C  0 

1  .6b 

3.46 

1  .52 

43 

275 

10 

0.0100 

1.68 

1715 

522 

10  700 

7  .3 

520 

11 

0  .0086 

1  .44 

2000 

101C 

39200 

2.7 

520 

7 

0 .0200 

3.36 

2960 

101C 

86000 

5  .82 

44 

275 

15 

0.00444 

0.75 

1095 

522 

43800 

2.98 

520 

14 

0 .0051 

0.86 

1540 

3  010 

23200 

1  .58 

520 

11 

0 .0083 

1  .39 

17  70 

1010 

30900 

2 .10 

45 

200 

172 

0.0000338 

0  .0057 

239 

12 

61 

18 

0.00309 

0.52 

513 

346 

22000 

1  .50 

122 

3  8 

0.00303 

0.52 

685 

490 

19600 

1  .33 

183 

16 

0 .0039 

0.65 

800 

6  0  C 

17800 

1.21 

305 

18 

0.00309 

0.52 

1025 

773 

17503 

1.19 

425 

16 

0 .00309 

0 .52 

800 

911 

7700 

0.52 

700 

16 

0.0039 

0.65 

1540 

1175 

17500 

1.19 

15 

61 

122 

183 

305 

425 

3  8 

38 

38 

38 

38 

0.00069 

0 .00069 

0  .00069 
0.00069 

0 .00069 

0.116 

0.116 

0 .110 
0.116 

0 .116 

217 

331 

376 

548 

640 

346 

490 

600 

773 

918 

3700 

4560 

3940 

5000 

49C0 

0.25 

0  .31 
0.27 

0  .34 

0.33 

19 

122 

48 

0 .000435 

0.073 

228 

490 

217  0 

0.15 

183 

46 

0.000475 

0.009 

388 

600 

4180 

0.28 

v2S 

39 

0.000658 

0.011 

662 

918 

5220 

0.35 

700 

39 

0.000658 

0.011 

800 

1175 

4650 

0 .32 

22 

6 1 

16 ,4 

0 .0037 

0.62 

386 

3  57 

1  1800 

o 

CO 

o 

122 

16.0 

0.0039 

0.65 

50  2 

490 

10600 

C  .72 

305 

13.8 

0 .00o22 

o 

• 

ro 

CO 

800 

700 

13000 

o 

CO 

xO 

425 

13.4 

0.00555 

0.93 

970 

918 

11200 

0.76 

112 


Table  8 


Comparison  of  Plastic  Flow  Stress  with  Control  Tests 


Series  No. 

Average 

Average 

0 

Con  - 

Tests 

6  Target 

Material 

Po 

( eq . 191 ) 

kg/cm2 

S 

max 

% 

(S  )* 

max 

P  s 

,  Fo  .  max 

( eq  .  1 94 ) 
kg/cm2  % 

( c 

' ^max ’ 

trol 

CBR 

% 

qu 

kg/cm 

31  Cambar  wax 

0 .43 

-23 

528 

1.62 

+  13 

169 

_ 

0  .511 

M-348 
(at  31 °C ) 

32  Estar  3 

0.023 

+  13 

169 

0  .061 

+  21 

441 

- 

0.045 

grease 
(at  2 8°C ) 

34  paraffin 

127.o 

-  1 

81  . 

2  00.0 

+  3 

9 

- 

134.0 

(at  2 9°C ) 

37  clay 

2.62 

±  0 

10.0 

-46 

2120 

2.1 

2.2 

38  clayey  sand 

2.11 

-  6 

36 

5  .9 

±  9 

81 

4.3 

3  .0 

silt 

39  beach  sand 

0  .089 

+  18 

324 

0  .30 

+  73 

5320 

- 

- 

41  Princeton 

3.41 

-23 

528 

14.4 

-42 

1760 

6.3 

4  .4 

clay 

42  sandy  clay 

0.85 

-33 

1085 

2.21 

+  57 

3240 

4  .0 

2  .  8 

<*3  crushed  rock 

2.13 

+  59 

3480 

5  .27 

+  49 

24  J 

18.5 

- 

44  fine  sand 

1  .00 

+  39 

1520 

2 .22 

+  3u 

1155 

5  .8 

- 

12  Princeton 

0.56 

+  16 

25  6 

1  .14 

-54 

2905 

2.8 

2 . 1 

red  clay 

15  Ottawa  sand 

0.12 

_  0 

0 

0 . 30 

-17 

289 

- 

- 

( dry ) 

19  beach  sand 

0  .009 

-26 

675 

0  .25 

-60 

3600 

- 

- 

w  =  l  .34  g/cm 

22  Princeton 

0.76 

+  22 

485 

0 .91 

-11 

121 

1  .9 

1  .4 

clay  (w«2.65% 

'l 

/ 

Zs2 

=  9167 

L  s2  - 

23610 

Zs2/n 

=  652 

Z*2' n  . 

1685 

mean  scatter 

:  sm 

=  25.5% 

Sm  = 

41% 

1 

Hardnes  number 

computed 

from 

the  penetr 

ation  number  N. 

2Failure  stress 

computed 

from  the  tensile 

stress 

at  2  3°C  . 

»  vpm  low  i  \  tyyj'  y  ■  •  * 

MW  »■  '■■»--  I  n  .  — m£t*t**^f*  . 
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by  themselves  and,  particularly  with  a  rather  extensive  test 
program,  they  must  be  expected  to  have  a  certain  scattering. 

This  is  particularly  so  if  one  considers  the  great  vari¬ 
ability  of  the  properties  of  natural  soils  from  point  to 
point.  Thus,  perhaps  a  major  part  of  the  scattering  might 
be  traceable  to  actual  local  variation  of  the  soil  properties 
in  situ. 

If,  for  example,  we  disregard  the  values  for  the.  scatter¬ 
ing  of  pQ  from  (191)  of  test  series  No.  43  (  i  .e  .  impact  on 
crushed  rock  which  very  much  depends  on  the  surface  configur¬ 
ation  at  the  exact  impact  point),  the  average  scattering  is 
reduced  immediately  to  20%.  For  a  material  such  as  crushed 
rock  especially,  it  is  surely  important  whether  the  impact 
energy  is  directly  conducted  aw.ay  by  a  rigid  grain  skeleton 
or  whether  the  impact  is  a  resilient  and  elastic  one. 

The  standard  control  test  like  CBR  and  unconfined  com- 
pression  test  were  not  carried  out  on  clean  sands.  The  termer 
is  very  difficult  to  do  because  the  exact  depth  of  penetration 
is  difficult  to  establish  and  the  latter  is  not  possible  at 
all  in  a  non-cohesive  material. 

For  the  target  materials  of  series  31,  32  and  34,  i  .e  . 
for  the  cambar  wax,  the  Estan  3  grease  and  the  paraffin,  the 
mechanical  properties  furnished  by  the  manufacturer  were 
accepted . 
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The  penetration  number  N  for  a  grease  or  wax  is  the  pene¬ 
tration  (in  tenths  of  a  millimeter)  that  occurs  after  five 
seconds  when  a  90°  cone  weighing  250  grams  is  applied  to  the 
grease  at  25°  centigrade. 

For  these  materials  one  can  thus  determine  a  hardness 
number  using  equation  (7)  giving: 

(a)  for  the  Cambar  wax 

P  a  hj ,)  o 

Ht  =  0.898  —Z  =  0.898  — - =  0.51  (kg/cm  ) 

L  d2  ( 0  . 7  )  2 

v.  b )  for  the  Estan  3  grease: 

Hr  =  0.898  — =  0.045  kg/cm2 

L  (2.35)" 

For  the  paraffin  (Essowax  5010)  the  m;  ’facturer  gives  a  ten- 

r\ 

sile  strength  (a:  23°C)  of  2 o  kg/cm  .  If  one  uses  this  value 

to  calculate  the  shear  strength  and  applies  the  corresponding 

value  to  the  Prandtl  bearing  capacity  problem,  one  obtains  a 

value  of  pQ  =  134  kg/cm  which  is  close  enough  to  the  expen- 

2 

mentally  found  value  of  about  200  kg/cm  particularly  if  one 
considers  that  the  calculated  value  is  based  on  a  static 
strength  and  test  whereas,  for  the  case  of  the  dynamic  loads 
during  the  penetrometer  impact,  the  observed  values  should  be 
higher  . 

In  general,  the  values  for  the  plastic  flow  pressure  po 
evaluated  with  equation  (191)  are  more  consistent  and  are  in 
surprisingly  good  agreement  witr.  the  standard  control  test 
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values,  i.e.  with  the  unconfined  compression  tests  or  the 
corresponding  equivalent  CBR  values.  The  pQ  values  evalu¬ 
ated  by  equation  ( 1 9 M )  give  in  almost  all  cases  two  to  three 
times  higher  values  .  The  author  has  not  been  able  to  find  an 
explanation  so  far  for  this  fact,  except,  as  was  mentioned 
before,  when  calculating  pQ  from  (194)  two  experimental  values 
have  to  be  introduced  which  may  lead  to  a  multiplication  of 
any  errors.  However,  the  values  from  (194)  are  at  least  in 
agreement  as  far  as  the  order  of  magnitude  for  pQ  is  con¬ 
cerned  and  they  are  also  more  or  less  consistent  within  any 
test  series. 


4 .  Wave  Propagation  Velocities 

The  instrumenting  of  the  target  soils  with  accel¬ 
erometers  nad  one  main  purpose:  to  measure  the  elastic 
properties  of  these  target  materials.  They  all  were  built 
in  at  such  a  distance  from  the  point  of  impact  that  in  their 
vicinity  the  resulting  stresses  and  deformations  would  remain 
within  the  elastic  range.  furthermore,  they  were  to  check 
the  validity  of  the  assumption  made  generally  that  the  energy 
carried  away  by  vibrations  may  be  neglected.  The  following 
values  were  evaluated  from  the  vibration  signals: 

vi 

(a)  for  the  propagation  velocity:  a^  -  — 
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(b)  for  the  period  of  vibration  T  =  —• 

(c)  maximum  amplitude  of 

Bi 

the  acceleration:  bi  max  = 

where  g  =  981  cm/sec.  is  the  gravitational  acceleration. 


The  next  page  demonstrates  the  evaluation  of  these  quan¬ 
tities  from  the  test  data.  The  values  obtained  are  listed  in 
Table  9.  From  the  average  values  for  "a"  the  corresponding  E 
moduli  were  found  according  to  equation  (145)  which  are  shown 
in  Table  10.  According  to  this  equation  (145)  they  can  be 
determined  by: 


E 


(1  +  v)(l  -  2v)  _  2 

f a 


(145a) 


=  j*  a4  (for  v  =  0) 

=  0 .83  j  a2  ( for  v  =  0.25) 

The  values  found  are  in  very  good  agreement  with  those 
given  by  Barkan.  Their  scattering  is  relatively  small. 

The  photographs  of  series  41  which  have  the  least  dis¬ 
tortion  and  interference  show  clearly  that  the  signals  in  the 
soil  have  the  form  given  by  equation  (178).  If  we  take  the 
displacement : 

„  -Bx  i  (ax  -  ut ) 
u  =  Be  e 

and  find  the  acceleration  by  double  differentiation  with 
respect  to  time,  one  gets: 
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Table  9 

Evaluation  of  the  Vibration  Signals 


Test 

Series 

S  page 

Drop 

Height 

H 

cm 

rl 

cm 

r2 

cm 

ai 

m/sec 

a2 

m/sec 

Period 

T1 

rn  *  s  0  c  • 

Period 

T 

2 

m . sec  . 

b, 

1  max 
cm/sec . ^ 

b. 

2  max 
cm/ sec . ^ 

37  -  A 

8 

275 

283 

482 

113 

105 

16  .3 

22.9 

26 .0 

8 . 8 

A 

9 

520 

290 

510 

151 

159 

13.5 

14  .0 

28.9 

11  .4 

38  -  A 

10 

275 

680 

1000 

400 

357 

14.6 

12.8 

63.2 

29  .0 

A 

11 

275 

680 

1000 

314 

323 

15.2 

13.5 

23.3 

8 . 8 

A 

12 

275 

680 

1000 

400 

454 

13.5 

13.0 

82.2 

36  .0 

A 

1  4 

520 

680 

1000 

272 

278 

13  .6 

14.6 

27.0 

1  "  .4 

39  -  A 

15 

275 

14G 

280 

200 

200 

12.0 

11.0 

38 .0 

A  /  .4 

A 

17 

520 

1  4  C 

280 

188 

200 

13.0 

12.0 

50  .8 

17.3 

A 

18 

520 

140 

280 

165 

170 

12.0 

6  .5 

108  .0 

28.9 

41  -  A 

19 

275 

244 

490 

244 

340 

7  .  / 

7 . 5 

40 .5 

15  9.0 

A 

20 

275 

242 

488 

242 

296 

5  .4 

7.0 

39  .0 

121.5 

A 

22 

520 

244 

491 

244 

230 

■ '  6.6 

8 .0 

57.8 

232.0 

42  -  A 

24 

275 

178 

396 

170 

208 

9  .6 

12.5 

57.8 

17.9 

43  -  A 

26 

275 

112 

233 

12'* 

123 

4  .5 

6  .1 

202  .0 

95  .4 

A 

27 

520 

107 

238 

153 

164 

6  .0 

6  .0 

179.0 

55.0 

44  -  A 

29 

275 

151 

375 

- 

417 

- 

10.0 

98  .1 

21.7 

A 

31 

520 

166 

394 

- 

329 

7  .5 

i  0  .0 

188.0 

41  .9 

r3 

a3 

T  3 

b3 

41  -  A 

19 

850 

283 

6  .1 

3.35 

A 

20 

825 

275 

7.25 

3.00 

A 

22 

630 

258 

6.25 

5  .00 

118 


x 

2 


) 

i 

4 

3 
9 
0 

5 
0 
9 

4 
0 
.7 
.9 


r\  •  ..  /  -  «  £  r>  \  • 

«•  _  £  XAVW  T  ip/  -i  COL 

u  =  B  to  e  e 


(195  ) 


i  .e  . 


except  for  the  factor  w 


2 


nothing  has  changed. 


Table  30 


Determinati 

on  of  E 

-modulus 

Test 
Series 
No  . 

Target 

Material 

P 

2  ,  4 

cm  /sec 

a 

ai/sec 

a2 

m2/sec2 

E ( v  =  0) 

kg/ cm^ 

E ( v  =  1/4) 

kg/cm^ 

37 

Clay  soil 

0 .00126 

135 

18,200 

230,000 

190 ,000 

38 

Clayey  sand 
-  silt 

0 .00138 

350 

122,000 

1  ,680  ,000 

1 ,390  ,000 

39 

Beach  sand 

0.00137 

188 

35,200 

480,000 

400  ,000 

41 

Princeton  clay 

0 .00156 

274 

45 ,000 

1 ,117,000 

970  ,000 

42 

Sandy  clay 

0.00141 

189 

36 ,000 

510  ,000 

422  ,000 

43 

Crushed  rock 

0  .0230 

141 

20 ,000 

460  ,000 

380  ,000 

44 

Fine  sand 

0.00153 

373 

139 ,000 

2,100  ,000 

1  ,740,000 

5  .  Evaluation  of  Vibration  Energy 


The  energy  furnished  by  a  force  P  moving  through  a 
distance  s  is  given  by: 


.  fpd../' 

Jo  Jo 


muds 


W 


P  ds 


(196) 
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If  we  determine  an  average  value  P  through  the  integral  inter¬ 
val  from  the  test  data,  then  the  total  vibration  energy  could 
be  : 


W  =  P 


av 


Considering  a  half-period  1/2  T  with  a  sinusoidal  change 
of  the  acceleration  from  zero  to  a  maximum  amplitude  b  and 
back  to  zero,  we  may  write  for  the  velocity  change: 


1/2  T 


2  .  i  b  i  T  =  i  £T 

•n  2  2  2  7T 


(197) 


Analogously  the  distance  traveled  during  the  half-period 


is  : 


„  2  1  1  bT‘ 

S~TIV2T  =  2 


and  the  average  force  becomes 


(19  8) 


P  =  m  ~ 
m  u 


where  m  designates  the  accelerated  mass. 

The  total  energy  thus  is  determined  by: 


b  2T  2 


W  = 


m 


2  TT 


(199) 


(  200  ) 


If  one  assumes  the  waves  propagate  spherically,  the  acceler¬ 
ated  mass  is : 


m  =  3 


it  r 


(  201  ) 
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and  the  total  energy  then  is: 


W 


r3b2T2 


{  202  > 


According  to  this  formulation  the  maximum  energy  of  vibration 
evaluated  for  all  test  series  was  found  for  series  38  (page 
A-10).  Using  the  corresponding  values  from  Tables  9  and  10, 
we  found  for  this  case: 


W, 


29  .5 


680  •  82.2  *  13.5' 


0.00138  =  18,500  cm.gm 


The  total  impact  energy  however  for  this  test  was: 


WQ  =  1573  •  275  =  433,000  cm.gm. 
giving  a  ratio  N  of  vibration  energy  to  impact  energy: 


v  _  18 ,500  _  4  3 

WQ  ”  433,000 


Since  this  was  by  far  the  largest  value  found,  most  other 
« 

values  being  a  significant  amount  below  this  ratio,  we  con¬ 
cluded  that  for  our  experiments  the  energy  dissipated  by  the 
vibrations  of  the  half-space  is  indeed  negligible  . 


r 
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X .  Summary  and  Conclusions 

The  purpose  of  this  investigation  was  to  explore  a  method 
that  would  permit  the  reliable  assessment  of  surface  soil 
properties  for  the  prediction  cf  deformations  during  and  after 
aircraft  operations.  After  the  historical  development  of  the 
penetrometer  concept  is  briefly  traced,  we  present  the  basic 
theoretical  aspects  of  the  problem.  A  penetrometer  with  a 
piezo-electric  accelerometer  as  the  main  feature  was  designed. 
This  penetrometer  was  tested  and  for  laboratory  experimenta¬ 
tion  purposes  evaluated.  While  the  instrument  as  designed 
certainly  should  not  yet  be  considered  an  ultimate  version  of 
a  new  operational  penetrometer,  the  principle  tested  and  the 
methods  employed  were  found  valid  and  useful.  The  modifica¬ 
tions  suggested  are  to  make  the  penetrometer  a  complete  sphere, 
use  three  mutually  perpendicular  accelerometers,  record  the 
signal  on  magnetic  tape  and  eliminate  the  cable  by  telemeter¬ 
ing  the  signal  through  a  built-in  transmitter  to  a  receiver, 
the  surface  of  the  sphere  acting  as  an  antenna.  With  these 
modifications  it  is  believed  that  a  system  will  be  achieved 
which  will  permit  the  remote  determination  of  the  mechanical 
strength  of  a  soil  surface  with  a  reliability  sufficient  for 
engineering  purposes. 

The  values  found  for  the  flow  pressure  were  in  reasonable 
agreement  with  the  control  test  values  and  were  consistent 


mmsz-  ,  ’ 11  *r 
«  ~ 

.  it 
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among  rhemselves.  The  solution  for  a  sphere 
four  parameter  half-space  was  developed  which 
a  rather  complicated,  non-linear  result  that 
evaluated  numerically.  Finally,  it  was  shown 
dissipated  by  vibrations  when  a  body  impacts 
for  the  velocities  and  material  investigated 
gible  . 


impacting  on  a 
gives,  however, 
must  yet  be 
that  the  energy 
on  a  half-space 
is  indeed  negli- 
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The  Aerial  Penetrometer 


(U.S.  Air  Force) 


(a) 


(b) 


Figure  2 

Cone  Penetrometer  of  U.S.  Army  Corps  of  Engineers 

(a)  field  instrument 

(b)  laboratory  instrument 
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Figure  8 

Deformations  in  the  Soil  after  CBR 
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Figure  il 

Ram  Hardness  versus  Compressive  Strength 


Figure  12 

Picture  of  Impact  Penetrometer 
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Schematic  lest  Setup 


Test 


Series 


No.  31 


Target  Material:  Cambar  Wax  M-348  at  31  C. 
Calibration  Factor  of  Accelerometer:  K  =  8.6  mV/g 


Drop  Height:  H  =  61  cm 
Oscilloscope  Calibration: 
M_  s  200  mV/cm 

a 

Mt  =  5  msek/cm 


t  -  0.16  cm 


■HBrnaiim 


Drop  Height*.  H  =  122  cm 
Oscilloscope  Calibration: 


Ma  *  200  raV/cm 


Mt  = 


5  msek/cm 


t  =  0.64  cm 
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Test  Series  No.  31 

Target  Material:  Cambar  Wax  M-3'i?  at  31°C. 
Calibration  Factor  of  Accelerometer:  K  =  8.6  mV/g 


Drop  Height:  H  =  183  cm 
Oscilloscope  Calibration: 
Ma  =  500  mV/cm 
M^  =  5  msek/cm 

t  s  0.8  cm 


Drop  Height:  H  =  305  cm 
Oscilloscope  Calibration: 
Ma  -  500  mV /cm 
=  5  msek/cm 

t  =  3 . 2  cm 
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MaHKZflC: 
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Series 


No.  31 


Target  Material:  Cambar  Wax  M-348  at  31°C. 
Calibration  Factor  of  Accelerometer:  K  =  8 


Drop  Height:  H  =  425  cm 
Oscilloscope  Calibration: 


Ma  =  500  mV/ cm 

M-t  =  5  msek/cm 

/ 

N 

t  =  3 . 8  cm 

4+44- 

4444 

444+ 

+H- 

\ 

J 

fflm 

A-4 


T1  ~  ~ 

A  O  U 


S  6  i‘  I  e  s  No.  32 


Target  Material:  Lubricating  Grease  Estan  3  at  28  C 
Calibration  Factor  of  Accelerometer:  K  =  8.6  mV/g 


Drop  Height:  H  =  61  cm 
Oscilloscope  Calibration: 


50  mV/cm 


Mt  =  10  msek/cm 


t  =  12.8  cm 


Drop  Height:  H  =  122  cm 
Oscilloscope  Calibration: 
Ma  =  50  mV/cm 

=  5  msek/cm 


t  =  15.0  c  rr. 


Calibration  Factor  of  Accelerometer:  K  =  8,6  mV/g 


Drop  Height:  H  =  183  cm 


Target  Material:  Paraffin  at  28°C. 
Calibration  Factor  of  Accelerometer:  K 


=  8.6  mV/g 


Drop  Height:  H  =  61  cm 
Oscilloscope  Calibration: 
Ma  =  2,000  mV/cm 


0.2  msek/cm 


0.1  cm 


Drop  Height:  H  =  122  cm 
Oscilloscope  Calibration: 
M  =  5,000  mV/cm 


Mt  =  0.2  msek/cm 


0.1  cm 


KlnHrS  I  IMlFI 


rr«"« 
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Test  Series  No. 37 

Target  Material:  Undisturbed  Clay  Soil 
Calibration  Factors:  Penetrometer:  K  =  8.6  mV/g 

r^ :  K  =  34  mV/g;  r 2 :  K  =  56  mV/g 


Drop  Height:  H  =  275  cm 
Oscilloscope  Calibration: 

Penetrometer : 

Ma  =  1000  mV/cm 
M+  =  10  msek/cm 


r^  =  283  cm 

Ma  =  0.5  mV/cm 

Mt  =  10  msek/cm 

r2  =  482  cm 

Ma  =  0.5  mV/cm 

Mt  =  10  msek/cm 
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Test  Series  No.  37 

Target  Material:  Undisturbed  Clay  Soil 
Calibration  Factors:  Penetrometer:  K  =  8.6  mV/g 

v±:  K  =  34  mV/g;  r2 :  K  =  56  mV/g 


Drop  Height:  H  =  520  cm 
Oscilloscope  Calibration: 


Penetrometer : 

Ma  =  1000  mV/cm 
M-t  =  10  msek/cm 

r^  =  290  cm 

Ma  =  0.5  mV/cm 

M-j-  =  10  msek/cm 


r2  =  510  cm 

M3  =  0.5  mV/cro 

Mt  =  10  msek/cm 

Penetrometer: 

Ma  =  2000  mV/cm 


Test 


A  - 1 1 


Series  No  38 

Target  Material:  Clayey,  Sandy  Silt,  Undisturbed 
Calibration  Factors:  Penetrometer:  K  =  8.G  mV/g 

rl :  K  =  31  mV/g;  r2 :  K  =  34  mV/g 


Drop  Height:  H  =  275  cm 
Oscilloscope  Calibration: 

Penetrometer : 

Ma  =  200  mV/cm 

M-(-  =  2  msek/cm 

=  680  cm 

Ma  =  1  mV/cm 

M-j-  =  10  msek/cm 

r2  =  1000  cm 
Ma  =  500  mV/cm 

Mt  =  10  msek/cm 


A-12 


Test  Series  No.  3 S 

Target  Material:  Clayey,  Sandy  Silt,  Undisturbed 
Calibration  Factors:  Penetrometer :  K  =  8.6  mV/g 

:  K  =  31  mV/g;  r^:  K  =  34  mV/g 


Drop  Height:  H  =  275  cm 


Oscilloscope  Calibration: 


Penetrometer : 


Ma  =  200  mV/cm 

Mt  =  10  msek/cm 


A 


i  ! 


i  I 


rn . 


r^  =  680  cm 

M  =  1  mV/cm 

Cl 

Ho.  =  10  msek/cm 


T2  -  1000  cm 
Ma  =  500  mV/cm 

M*.  =  10  msek/cm 


Test 


Series 


No.  38 


Target  Material:  Clayey,  Sandy  Silt,  Undisturbed 
Calibration  Factors:  Penetrometer:  K  =  8.6  mV/g 

r±:  K  =  31  mV/g;  r2 :  K  =  34  mV/g 


M-j-  “  200  rasek/cm 


A-14 


Test  Series  No.  3B 

Target  Material:  Clayey,  Sandy  Silt,  Undisturbed 
Calibration  Factors:  Penetrometer:  K  =  8.6  mV/g 

r^ :  K  =  31  mV/ g;  r2 :  K  =  34  mV/g 


Drop  Height:  H  =  520  cm 


Oscilloscope  Calibration: 

Penetrometer : 

Ma  =  500  mV/cm 

Mt  =  2  msek/cm 

r^  =  680  cm 

Ma  =  500  mV/cm 

Mt  =  20  msek/cm 


r2  =  1000  cm 
Ma  =  200  mV/cm 

Mt  =  20  msek/cm 


tawacassL.  5»y  a— i 

"  ''r.V'*Shz?" J  ' 


Test 


Series 


No.  3Q 


Target  Material:  Beach  Sand,  Freely  Deposited 
Calibration  Factors:  Penetrometer:  K  =  8.6  mV/g 

r1 :  K  =  31  mV/g;  r K  =  34  mV/g 


Drop  Height:  H  =  275  cm 


Oscilloscope  Calibration 


Penetrometer : 


Ma  =  100  mV /cm 


=  5  msek/cm 


r^  =  140  cm 
M  _  =  500  mV/cm 
Mt  =  10  msek/cm 


V2  -  280  cm 
Ma  =  500  mV/ cm 
Mt  =  10  msek/cm 


■VtilMtt&H 

SlViUWiififij 


A- 


Test  Series  No. 39 

Target  Material:  Beach  Sand,  Freely  Deposited 
Calibration  Factors:  Penetrometer:  K  =  8.6  mV/g 

rl!  K  =  31  mV/g;  r2:  K  =  34  mV/g 


Drop  Height:  H  =  520  cm 


Oscilloscope  Calibration: 

Penetrometer : 

Ma  =  200  mV/cm 
=  5  msek/cm 

r^  =  140  cm 
Ma  =  2  mV/cm 

=  10  msek/cm 


r2  =  280  cm 
Ma  =  1  mV /cm 

Mt  =  io  msek/cm 


SR55 


Target  Material:  Princeton  Red  Clay 
Calibration  Factors:  Penetrometer:  K  =  8.6  mV/g 

rl :  K  =  34  mV/g;  r2 :  K  =  34  mV/g 

r3 :  K  =  170  mV/g 

Drop  Height:  H  =  275  cm 
sciiloscope  Calibration* 


Penetrometer : 
Ma  =  500  mV /cm 


M-t  ~  5  msek/cm 


=  244  cm 
Ma  s  500  mV/cm 
=  10  msek/cm 


r2  =  490  cm 
Ma  =  1  mV/cm 

Mt  =  10  msek/cm 


r„  =  850  cm 


=  200  mV/cm 
=  5  msek/cm 


Teat  S  e  r  i  e 


s  No.  41 


Target  Material:  Princeton  Red  Clay 
Calibration  Factors:  Penetrometer:  K  =  8.S  mV/g 

r1 :  K  a  34  mV/g;  r2 :  K  a  34  rov/g 

r3:  K  =  170  mV/g 


Drop  Height:  H  =  275  cm 
Oscilloscope  Calibration: 


Penetrometer : 


M,  =  500  mV/cm 


5  msek/cm 


=  242  cm 
a  500  mV/cm 
4*  =  10  msek/cm 


^2  -  488  cm 
4a  =  1  mV/cm 

a  10  msek/cm 


rg  =  825  cm 
Ma  a  0.2  mV/cm 
Mt  *  5  msek/cm 


Target  Material:  Princeton  Red  Clay 
Calibration  Factors:  Penetrometer:  K 

r,  :  K  -  34  mV/g; 
r3:  K  =  170  mV/g 

Drop  Height:  H  =  275  cm 
Oscilloscope  Calibration: 


T  e  8  t 


Series 


No.  41 


Target  Material:  Princeton  Red  Clay 
Calibration  Factors:  Penetrometer:  K  =  8.6  mV/g 

r^:  K  =  34  mV/g;  :  K  =  34  mV/g 

r3:  K  =  170  mV/g 

Drop  Height:  H  -  520  cm 
Oscilloscope  Calibration: 

Penetrometer : 

=  500  mV/cm 

d 

Mt  =  5  msek/cm 


r.  =  244  cm 
1 

Ma  =  1  mV/ cm 

Mt  =  10  msek/cm 


r0  =  491  cm 
Ma  =  2  mV/cm 

=  10  msek/cm 

r3  =  830  cm 
Ma  =  0.2  mV /cm 
M*.  =  5  msek/cm 


Target  Material:  Silty 
Calibration  Factors:  Pe 


rl: 

Drop  Height:  H  =  275  cm 
Oscilloscope  Calibration: 

Penetrometer : 

M  a  =  200  mV/cm 
M  •(-  =  2  msek/cm 

r^  =  198  cm 
Ma  =  1  mV/cm 

Mt  =  10  msek/cm 

r2  =  415  cm 
M-  =  500  mV/cm 


10  msek/cm 


(Plowed  Field) 


etrometer:  K  =  8.6  mV/g 


K  a  34  mV/g;  r2 : 


K  =  34  mV/g 


lest 


b  e  r  i  e  s 


Target  Material:  Silty  Sand  (Plowed  Field) 

Calibration  Factors:  Penetrometer:  K  =  8.6  mV/g 

rjL:  K  =  34  mV/g;  v ^ :  K  =  34  mV/g 


Drop  Height:  H  =  275  cm 


Oscilloscope  Calibration: 


Penetrometer : 

Ma  =  200  mV/cm 
=  2  msek/cm 

r^  =  178  cm 
Ma  =  500  mV/cm 
Mt  =  10  msek/cm 


r2  =  396  cm 
Ma  =  200  mV/cm 
M-f  =  10  msek/cm 


i  c  o  t  o  c  r  i  c  5  No*  4  2 


Target  Material:  Silty  Sand  (Flowed  Field)  ' 

Calibration  Factors:  Penetrometer:  K  =  8.6  roV/g 

rl :  K  =  34  mV/g;  r2 :  K  =  34  mV/g 


Drop  Height:  H  =  520  cm 
Oscilloscope  Calibration: 

Penetrometer : 

Ma  =  200  mV/cm 
M^.  =  2  msek/cm 


r^  =  178  cm 
Ka  =  1  mV/cm 

M^.  =  10  msek/cm 


*2  =  396  cm 
Ma  =  500  mV/cm 


M-j-  s  10  msek/cm 


Test 


Series 


No.  *‘l  3 


T  e  s 


Seri 


Target  Material:  Crushed  Rock,  Compacted 
Calibration  Factors:  Penetrometer:  K  =  8.6  mV/g 

r±:  K  =  34  mV/g;  r2 :  X  =  34  mV/g 


Drop  Height:  H  =  520  cm 

Oscilloscope  Calibration 

Penetrometer : 

Ma  =  500  mV/cm 
M.J.  =  1  msek/cm 


r1  -  107  cm 
Ma  =  2  mV /cm 

M.J.  =  10  msek/cm 


F  'W'‘ : 


• ■ 


r2  '  =  238  cm 
=  1  :  V/cm 

(4*.  =  10  msek/cm 


Test  Series  No,  43 

Target  Material:  Crushed  Rock,  Compacted 

Calibration  Factors:  Penetrometer:  K  =  8.6  tnV/g 

r1:  K  =  34  mV/g;  r2 :  K  -  34  mV/g 

Crop  Height:  H  =  520  cm 
Oscilloscope  Calibration: 

Penetrometer : 

Ma  =  1000  mV/cm 
Mt  =  1  msek/cm 


r^  a  107  cm 
Ma  =  2  mV/cm 

Mt  a  10  msek/cm 

r2  =  238  cm 

K  =  1  mV /cm 

M+.  =  10  msek/cm 


Test 


Series 


No.  44 


Target  Material:  Uniform,  Fine  Sand,  Natural  Deposit 
Calibration  Factors:  Penetrometer:  K  =  8.6  mV/g 

r1:  K  =  34  mV/g;  r2 :  K  =  34  mV/g 


Drop  Height:  H  =  275  cm 
Oscilloscope  Calibration: 

\ 


Penetrometer : 

Ma  =  200  mV/cm 

M t  =  2  msek/cm  j 

i 

t 

r^  =  151  cm 
Ma  =  1  mV/cm 

Mt  =  10  msek/cm 


r2  =  375  cm 
Ma  =  500  mV/cro 
Mt  =  10  msek/cm 


A-3C 


Test  Series  Ho  4h 

Target  Material:  Uniform*  Fine  Sand*  Natural  Deposit 
Calibration  Factors:  Penetrometer:  K  =  8.6  mV/g 

r^ :  K  =  34  mV/g;  •  K  =  34  mV/g 

Drop  Height:  H  =  520  cm 
Oscilloscope  Calibration: 

i 

Penetrometer: 

Ma  ~  500  mV/cm 
Mt  =  2  msek/cm 

r.  =  160  cm 
1 

Ma  =  2  mV/cra 

=  10  msek/cm 


rj  =  337  cm 
-  5  mV/cm 

Q 


j  '  ^  » v  ' 


10  msek/cm 


-30 


A- 31 


c  « 


2  N  o «  4  4 


Target  Material:  Uniform,  Fine  Sand,  Natural  Deposit 
Calibration  Factors:  Penetrometer:  K  =  8.6  mV/g 

r^:  K  =  34  mV/g;  r2  :  K  =  34  mV/g 


Drop  Height:  K  =  520  cm 
Oscilloscope  Calibration: 

Penetrometer : 


Mt  -  10  msek/cm 


r2  =  394  cm 
Ma  =  500  mV/cm 
=  10  msek/cm 


7 


Oi 


—Ai  n * 


r 


Test  Series  No.  45 
Target  Material:  Swamp 

Calibration  Factors:  Penetrometer:  K  =  8.6  mV/g 

rl :  K  =  34  mV/g;  r2 :  K  =  34  mV/g 


Drop  Height:  H  =  275  cm 


Oscilloscope  Calibration: 

/ 


Penetrometer : 

=  50  mV/cm 

Mt  =  20  msek/cm 


r^  =  76  cm 

Ma  =  500  mV/cm 
Mt  =  20  msek/cm 


r2  =  76  cm 

Ma  =  500  mV/cm 
=  20  msek/cm 


I 

I 

I 

I 


Test 


Series 


No.  12 


V 


Target  Material:  Sandy  Silt  with  some  Clay; 

Natural  Surface  Deposit 

Calibration  Factor  of  Accelerometer:  K  =  8.6  mV/g 


.  rop 

He 

ight : 

H  =  183  cm 

Osci 

llo 

scope 

Calibration : 

= 

20 

o 

3 

< 

cm 

Mt  = 

2  mse 

k/cm 

t  = 

3  . 

7  cm 

Drop  Height:  H  =  305  cm 
Oscilloscope  Calibration: 
Ma  =  200  mV/cm 
Mt  -  2  msek/cm 

t  =  4 . 5  cm 


-2 


Test  Series  No. 12 

Target  Material:  Sandy  Silt  with  some  Clay; 

Natural  Surface  Deposit 

Calibration  Factor  of  Accelerometer:  K  ~  8.6  mV/g 

Drop  Height:  H  =  425  cm 
Oscilloscope  Calibration: 

Ma  =  500  mV/cm 
=  2  msek/cm 

t  =  7 . 9  cm 


Drop  Height:  H  =  700  cm 
Osc  lloscope  Calibration: 
M_  =  500  mV/cm 

Q 

Mt  =  5  msek/cm 

t  ~  5  -  2  cm 


o 


\T~ 
11W  « 


A  C 
J.  -» 


fP  A  ^ 

i  ^  o  v.  c  c;  i  x  c  o 


Target  Material:  Ottawa  Sand,  U.S.  Sieve  Size  No.  30-40 
Calibration  Factor  of  Accelei ometer :  K  =  8.6  mV/g 


Drop  Height:  H  =  183  cm 
Oscilloscope  Calibration: 
Ma  =  100  mV/cra 
Mt  =  5  msfck/cm 

t  =  5 . 7  cm 


Drop  Height:  H  =  305  cm 
Oscilloscope  Calibration: 
Ma  =  200  mV/cm 

Q 

M^  =  5  msek/cm 


t 


6.4  cm 


T  «  8  t  Series  Ho^lS 

Target  Material:  Ottawa  Sand,  U.S.  Sieve  Size  No.  30-40 
Calibration  Factor  of  Accelerometer:  K  =  8.6  mV/g 


Drop  Height:  H  =  425  cm 
Oscilloscope  Calibration: 
Ma  =  200  mV/cm 
Mt  a  5  msek/cm 


t 


6.8  cm 


Target  Material:  Beach  Sand,  Natural  Deposit 
Calibration  Factor  of  Accelerometer:  K  =  8.6  mV/g 

Drop  Height:  H  =  61  cm 
Oscilloscope  Calibration: 

Ka  =  50  mV/cm 

Mt  =  5  msek/cm 

t  =  6 . 8  cm 

i. 


Di'op  Height:  H  -  122  cm 
Oscilloscope  Calibration: 
Ma  =  50  mV/cm 

Mt  =  5  msek/cm 

t  =  7 . 0  cm 


Target  Material:  Beach  Sand,  Natural  Deposit 
Calibi'ati 'n  Factor  of  Accelerometer:  K  =  8 , 6 


Drop  Height:  H  =  183  cm 
Oscilloscope  Calibration: 
M_  =  100  mV/cm 

a 

=  10  msek/cm 

t  =  7 . 0  cm 


Drop  Height:  H  =  30,6  cm 
Oscilloscope  Calibration: 
Ma  =  200  mV/cm 
Mt  =  5  msek/cm 


mV/g 


t 


7,5  cm 


Target  Material:  Beach  Sand,  Natural  Deposit 
Calibration  Factor  of  Accelerometer:  K  =  8.6  mV/g 


Drop  Height:  H  -  425  cm 
Oscilloscope  Calibration: 
M.  =  200  mV/cm 

a. 

Mt  =  5  msek/cm 

t  =  7 . 8  cm 


\ 


Drop  Height:  H  =  700  cm 
Oscilloscope  Calibration: 
M_  =■  2  00  mV /cm 

a 

Mt  =  5  msek/cm 


4* 


t 


8.0  cm 


Target  Material:  Princeton  Red  Clay; 

Moisture  Content  ~  2.65%,  w  =  86.8  lb/ft' 

Calibration  Factor  of  Accelerometer:  K  =  8.6  mV/g 


Drop  Height:  H  =  61  cm 
Oscilloscope  Calibration: 


Ma  =  100  mV/cm 
Mt  -  2  msek/cm 

t  =  2 . 3  cm 


4  I 

_* 


Drop  Height:  H  =  122  cm 
Oscilloscope  Calibration: 
M_  =  100  mV /cm 

a 

Mt  =  2  msek/cm 


t  =  3 . 0  cm 


7 


Target  Material:  Princeton  Red  Clay; 

Moisture  Content  =  2*65%,  w  =  86*8  lb/ft 

Calibration  Factor  of  Accelerometer:  K  =  8.6  mV/g 


Drop  Height:  H  =  183  cm 
Oscilloscope  Calibration: 


Ma  =  200  mV/cm 


Mt  =  2  msek/cm 


t  =  4 . 0  cm 


I 


Drop  Height:  H  =  305  cm 
Oscilloscope  Calibration: 


Ma  =  200  mV/cm 


M+  =  2  msek/cm 


t  ~  4.2  cm 


Test  Series  No.  22 

Target  Material:  Princeton  Red  Clay; 

Moisture  Content  -  2.65%,  w  =  86.8  lb/ft 

Calibration  Factor  of  Accelerometer:  K  =  8.6  mV/g 

Drop  Height:  H  =  425  cm 
Oscilloscope  Calibration: 

Ma  =  500  mV/cm 
Mt  s  2  msek/crn 
t  =  4.8  cm 


Drop  Height:  H  =  122  cm 
Oscilloscope  Calibration: 

Ha  =  200  mV/cm 
Mt  =  2  msek/cm 

t  -  3.0  cm  f - 1  r  i  I  T'“T 
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13.  abstract  The  report  presents  part  of  the  scientific  investigations  carried  out 
under  a  task  program  related  to  the  determination  and  assessment  of  surface 
dynamic  soil  properties  by  laboratory  instrumentation.  It  first  reviews  the  his¬ 
torical  development  of  the  hardness-penetrometer  concept  and  presents  a  pilot  model 
for  a  new  aerial  penetrometer.  Basic  theoretical  considerations  are  then  advanced 
regarding  the  stress-strain  relationships  in  axial  symmetry  and  the  dynamic  loads 
associated  with  the  drop-impact  penetrometer.  The  existing  theory  of  elasticity 
solutions  are  then  reviewed  and  the  solutions  for  a  sphere  impacting  on  an  ideally 
plastic  and  on  an  elasto-plastic  semi-infinite  solid  are  presented.  The  Hertz 
solution  of  an  elastic  sphere  impacting  on  an  elastic  half-space  is  modified  by 
using  the  elastic-viscoelastic  analogy  to  find  a  solution  for  a  sphere  dropping  onto 
a  four  parameter  viscoelastic  half  space. 

The  problem  of  energy  transfer  during  impact  by  vibrations  is  discussed  shortly 
The  data  for  15  test  series  carried  out  with  the  pilot  model  penetrometer  on  natural 
soils  and  on  some  artifical  materials  are  presented  and  evaluated.  It  is  concluded 
that  the  principle  of  the  drop  impact-penetrometer  is  applicable  in  experiments 
for  assessing  the  mechanical  properties  of  a  natural  soil  surface.  The  proposed 
modification  of  a  prototype  penetrometer  should  also  increase  its  capability  and 
reliability  as  a  remote  sensor  for  operational  applications.  The  instrument  as 
envisaged  will  not  only  permit  determination  of  mechanical  surface  properties  on 
terrestrial  soils  but  should  also  be  extremely  useful  for  the  exploration,  by 
remoted  techniques,  of  lunar  and  other  extraterrestrial  surfaces. 
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